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We study theoretically the return probability experiment, used to measure the dephasing time Tj, in a double 
quantum dot (DQD) in semiconducting carbon nanotubes (CNTs) with spin-orbit coupling and disorder induced 
valley mixing. Dephasing is due to hyperfine interaction with the spins of the '^^C nuclei. Due to the valley 
and spin degrees of freedom four bounded states exist for any given longitudinal mode in the quantum dot. 
At zero magnetic field the spin-orbit coupling and the valley mixing split those four states into two Kramers 
doublets. The valley mixing term for a given dot is determined by the intra-dot disorder and therefore the 
states in the Kramers doublets belonging to different dots are different. We show how nonzero single-particle 
interdot tunneling amplitudes between states belonging to different doublets give rise to new avoided crossings, 
as a function of detuning, in the relevant two particle spectrum, crossing over from the two electrons in one 
dot states configuration, (0,2), to the one electron in each dot configuration, (1, 1). In contrast to the clean 
system, multiple Landau-Zener processes affect the separation stage and the joining stage of each single-shot 
measurement and, even in a spin-orbit dominated situation, they affect the outcome of the measurement in a 
way that strongly depends on the initial state. We find that a well-defined return probability experiment is 
realized when, at each single-shot cycle, the (0,2) ground state is prepared. In this case, valley mixing increases 
the saturation value of the measured return probability, whereas the probability to return to the (0,2) ground 
state remains unchanged. Finally, we study the effect of the valley mixing in the high magnetic field limit; 
for a parallel magnetic field the predictions coincide with a clean nanotube, while the disorder effect is always 
relevant with a magnetic field perpendicular to the nanotube axis. 
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I. INTRODUCTION 

Since the Loss and DiVincenzo proposal,''^ electrons con- 
fined in quantum dots have become one of the most attractive 
platforms for realizing qubits in condensed matter systems-l^l 
For dots devised in GaAs based 2-dimensional electron gases 
(2DEGs), the spins of the confined electrons are affected by 
the hyperfine interaction with the nuclear spins — the nuclear 
spin /() -3/2 is common to the *^Ga, ^'Ga, and ''^As isotopes. 
In the pursuit of fault-tolerant quantum computation,^ this in- 
teraction has been recognized as one of the primary sources 
of dephasing,!^ and the challenge of avoiding this problem 
has led to the implementation and design of a number of ap- 
proaches such as (echo techniques, dynamic nuclear polariza- 
tion, etc.) with general success.^ ^ 

In parallel, there have been an increasing interest in devis- 
ing quantum dots in systems that can be isotopically pur ified 
leaving only spinless isotopes. This is the case of silicon,EfiE] 
graphene, and carbon nanotubes (CNTs), which is the sys- 
tem studied in this work. As a carbon based system, CNTs 
profit from the 99% natural abundance of the zero nuclear spin 
isotope, '^C. In particular, the gap found in semiconducting 
CNTs allows for confining electrons using gate-defined quan- 
tum dots; this has motivated a great deal of experimentall^^ES 
and theoretical studies.l22lEI]in particular, the generality of the 
presence of the spin-orbit coupling effects in quantum dots, 
due to the nanotube curvature,''^ has only recently been ex- 
perimentally recognized,'^ the effect persists even in disor- 
dered quantum dots occupied by hundreds of electronsj '^ l '^ l 
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FIG. 1. Double quantum dot in a CNT and the return probability ex- 
periment, (a) Scheme of the DQD, gate voltages Vl and Vr allow for 
the control of the number of electrons in each dot. (b) Single-shot ex- 
perimental cycle in gate voltage space consisting in preparing at point 
"p", evolving at point "s" and measuring at point "m". The detun- 
ing, e, is zero at the frontier between the (1,1) and (0, 2) regions, (c) 
Simplification of the applied pulse of detuning as a function of time 
during a single-shot experiment, (d) Sketch of the return probability, 
P(Tj), obtained from averaging many single-shot measurements, the 
saturation value, P„, is qualitatively presented. 

One of the experimental tools for measuring the dephasing 
is the return probability experiment (RPE), i.e., the measure- 
ment of the dephasing time, see FigjT] which involves both 
control and read-out of the charge state of a double quantum 
dot, (Nl, Nr), where Nl and Nr are the number of electrons in 
the left and the right dot, respectively. As sketched inside the 
charge stability map in gate-voltage space in FigjTJb), each 
single-shot measurement consists of five stages: (i) prepara- 
tion, the electron is prepared at the (0, 2) region; (ii) separa- 
tion, the system is taken adiabatically to the (1,1) region; (iii) 
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evolution, electrons are left to evolve under different environ- 
ments during a time Tj; (iv) joining, the system is taken adi- 
abatically back to the (0,2); and (v) measuring, the outcome 
is nonzero and set to 1 only if the measured charge state has 
returned to be (0, 2). For each t, many single-shot measure- 
ments are repeated and their outcomes are averaged leading 
to the return probability shape, Pijs), that characterizes the 
dephasing in the system, with a typical P(tj) shape shown 
in FigjTJd). The time scale for performing these fives stages 
should be much shorter than the inelastic relaxation time, T\ . 

In order to understand the impact of the hyperfine interac- 
tion, quantum dots have been devised in nanotubes with natu- 
ral abundance as well as in samples grown from 99% '^C (i.e., 
nuclear spin 1/2) enriched methaneP^l :„ Ref. [Tithe RPE 
was used to measure in a sample rich in '^C at zero mag- 
netic field, but nor finite magnetic field measurements neither 
results for samples rich in '^C are currently available. 

In absence defects and impurities, except from the hy- 
perfine interaction effects, there is no scattering between 
states from inequivalent valleys of the graphene-based band- 
structure. Already in such a clean limit the RPE in a nanotube 
double-dot leads to a rich variety of scenarios. Due to the spin 
and valley degrees of freedom a CNT DQD can be prepared in 
six different (0, 2) states and, once separated, in the evolution 
stage, the system have sixteen (1,1) states available for de- 
phasing. The situation is different from GaAs based DQDs in 
which case only the spin singlet state can be prepared and four 
(1,1) states are available at the evolution stage.E^MHin Ref.l44l 
for the RPE in clean nanotubes, we found nine dynamically 
inequivalent situations and five different lowest bounds for the 
saturation value of the return probability. The PRE outcome is 
highly dependent on the prepared state, the curvature induced 
spin-orbit coupling, and the diamagnetic and the Zeeman ef- 
fects of an external magnetic field. However, the lowest pos- 
sible saturation return probability. Poo (see Fig[TJd)), is 1/3, 
as in spin-only DQDs systems, which is above the 0.17 re- 
ported in the experiment of Ref. [14] Therefore, in order to 
expand our understanding of the RPE to non-ideal samples, 
we present here a new study that includes a spin-conserving 
valley mixing interaction which is induced by non-magnetic 
impurities and defects in the nanotube. 

At zero magnetic field, the valley mixing, which depends on 
the disorder profile and therefore is different for the two quan- 
tum dots, demands us to adopt a more general treatment than 
the clean case.l^ This arises because the products of triplet 
and singlet functions in valley and spin spaces are no longer 
the eigenstates of the (0,2) and the (1, 1) charge states. Each 
single dot eigenfunction is a particular combination of the two 
valleys and the single-particle tunneling amplitudes between 
solutions of different dots become non-trivial. Therefore, no 
simple selection rules apply to the mixing between the (0, 2) 
eigenstates and the (1,1) eigenstates, and multiple avoided 
crossings appear in the two-particle energy versus detuning. 

Moreover, we show that multiple Landau-Zener (LZ) pro- 
cesses are involved at the separation and the joining stages of 
the RPE and that the outcome of the experiment is highly de- 
pendent on the shapes of the detuning pulses at those stages — 
see the arrow starting at point "p" ("e") and ending at point "e" 



("m") in Figjlp^) for the separation (joining) stage and the as- 
sociated shape of the detuning pulse in Fig{T|c). When chang- 
ing the detuning the passage through multiple Landau-Zener 
processes is avoided only if the (0, 2) ground state is prepared. 
The return probability experiment becomes well defined be- 
cause the preparation guarantees that after the separation stage 
(at the beginning of the evolution stage) the electrons are in a 
(1, 1) state, i.e., they are really separated. In this situation, 
we find that the saturation return probability. Poo, is increased 
above 1/3; the hyperfine interaction behaves qualitatively as 
in the clean case whereas the (1,1) triplet-like states have an 
increased probability to return to (0, 2) due to a valley-mixing 
induced avoided crossing that affects the joining stage. 

We note that if the highest excited (0, 2) state is prepared, 
Poo ~ 1/6 is found for a broad region in parameters space. 
Surprisingly, this number is in good agreement with the satu- 
ration value reported in the experiment of Ref. [14] However, 
in our calculations, the multiple LZ processes also affect the 
separation and joining stages, generating a return probability 
less than one (even if no dephasing is allowed) by returning 
to the measurement point without waiting at the (1,1) region, 
and we find Pq = P(ts = 0) ~ 1/2. Despite the fact that the 
measurements at short times, P{Ts — » 0)~ 1, are the ones with 
greatest error bar we conclude that valley mixing does not ex- 
plain the reported experimental result. In order to test further 
the multiple Landau-Zener model, more measurements con- 
trolling both the prepared state — unknown in Ref. J4j — and 
the pulse shapes would be desirable. 

The behavior for nonzero magnetic fields is also unexplored 
experimentally. Here we show that in the limit of high parallel 
magnetic field the presence of valley mixing leads to the same 
predictions obtained for the clean case. This occurs because 
the diamagnetic effect dominates, so that the single-dot single- 
particle solutions are valley polarized and the mixing between 
valleys induced by disorder becomes negligible. On the other 
hand, a perpendicular external magnetic field, irrespective of 
its strength, is unable to avoid the effect of the valley mixing 
because it only introduces Zeeman interaction. Here again, 
the return probability experiment is not well defined due to 
the presence of multiple LZ processes affecting the separation 
and joining stages. We point out that this highly disorder- 
dependent mixing between the (1,1) and the (0, 2) states in a 
perpendicular magnetic field can be relevant for understand- 
ing Pauli Blockade measurements. 

The paper is organized as follows, in the next section we 
study the model for the single dot including the effective hy- 
perfine field interaction for a Kramers doublet and we intro- 
duce the two particle (0,2) and (1, 1) eigenstates. We show 
the effect of interdot tunneling and classify the resulting two- 
particle spectrum in Section III. Our return probability calcu- 
lation method and the results are presented in Section IV with 
conclusions in Section V. 



II. ISOLATED QUANTUM DOTS 

We include Coulomb effects in the constant interaction 
model, which is valid when the size quantization energy in the 
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quantum dot is large (short dots ) or w hen a strong dielectric 
substrate screens the interactionpSEH and Wigner molecular 
states and other interaction effects are not considered. Exper- 
imentally the absence of Coulomb exchange effects has been 
corroborated in many studies.'^ Our model is then con- 

structed from two particle Slater determinants based, for con- 
venience, on the single dot eigenfunctions of the left and the 
right dot. In this section we focus on those single dot wave- 
functions, on the effective hyperfine interaction seen by the 
solutions in a single dot, and we present the two-particle solu- 
tions assuming isolated dots. 



A. Dot Hamiltonian in valley and spin spaces 

The gate defined quantum dots we study are devised in 
semiconducting tubes with a bandgap that is due to e ither chi- 
rality or, for nominally metallic tubes, to curvature .'^SEl xhe 
low energy 7r-band electrons in the CNT are described by 
two gapped Dirac equations in one spacial dimensions (two 
Hamiltonian operators acting on 4-dimensional spinors due to 
the spin and pseudo-spin [weight in the two inequivalent sub- 
lattices] degrees of freedom), one for each of the two inequiv- 
alent valleys in reciprocal space, K and K' . Each of the two 
Dirac equations include the Zeeman interaction, the diamag- 
netic effect, and also spin-orbit coupling terms. For the quan- 
tum dots in nanotubes, since the electrostatic confinement po- 
tential that defines each dot is smooth on the length scale of 
the graphene honeycomb lattice spacing, the measured val- 
ley mixing effects are not likely to be generated by roughness 
at interfaces as opposite to other higher dimensional and/or 
etched-defined systems.'^Instead, as shown in Ref. 45 , valley 
mixing in nanotubes is expected as a result of the existence of 
nonmagnetic impurities and defects within each dot region. 

We work with an effective Hamiltonian for the lowest en- 
ergy bounded electron state. Such a description follows af- 
ter solving the Dirac equation,^^ and has been shown to fit 
well experimental results.'^ In contrast to quantum dots in 
2DEGs here, due to the coexistence of spin and valley degrees 
of freedom, the description of a bounded state in the dot is 4- 
dimensional. We introduce the identity and Pauli matrices in 
spin space cr„ with i - 0, x, y, z and the three dimensional spin 
vector a = (cr^, <Ty, cr-). The spin projection along the tube 
axis (z-direction) is denoted by cr -1, i (or alternatively its 
numerical version cr = ±1). Similarly, the identity and Pauli 
matrices in valley space are tj, with j = 0, 1,2,3, where we 
choose T = K, K' as the positive and negative projections of T3, 
respectively; then the 3-dimensional valley vector operator is 
just T = (ti,T2,T3). Neglecting the hyperfine interaction, the 
description for the right (R) and the left (L) quantum dots is 
given by the Hamiltonians: 

l^gs (B ■ cr) To -H goi-b (B ■ z) 0-0 raj , (1) 

with ^ = L, R. The first term describes the effect of the gate 
voltage applied to the dot ^, we assume that it only introduces 



a global energy shift of energy e^. The second term, is the 
spin-orbit coupling splitting energy A^^ between Kcr and K'cr 
states, where cf stands for the spin projection opposite to cr. 
Note that the curvature does not depend on the position in 
the nanotube and therefore the spin-orbit splitting seen by the 
lowest lying estate is the same for the two dots, A^„ = Af^,, if 
the dots have the same length.'^ The third term describes the 
valley mixing in the dot ^. For convenience we work with a 
single energy parameter (and its phase) to quantify the valley 



mixmg: 



A^ = A^ 



•iA' 



KK',2\ 



, = Arg{A' 



^ - iA^ ! 

KK'.\ 



(2) 



Finally, the fourth term describes the effect of an external 
magnetic field B, it includes both the usual Zeeman interac- 
tion (oc gs) and the diamagnetic effect (oc gojj,) which only 
appears if the external magnetic field has a non-zero compo- 
nent along the tube's axis. Typically the orbital gyromagnetic 
factor ^oib is greater than the spin gyromagnetic factor lead- 
ing to anisotropic magnetic effects.l^ '^l^^^ lft is useful to define 
the Zeeman and orbital splitting energies as Ef^ = gsi^glRl, and 
^orb = 2gorbjUB (B ■ z), respectively. 

For zero magnetic field (and for magnetic fields along the 
tube axis) the spin projection cr is a good quantum number, 
then we decompose the Hamiltonian of Eq. ([T]i as. 



where. 



T, 



Eg- = + o-Es/2 , 

= 7F i^KK' COS fiic' ' 4if ' sin V^K' ' ^o) , 

6i = cr Ai, - £orb , = ^J{A^^^,f + {sif 



(3) 

(4a) 
(4b) 

(4c) 



Note that for a given cr, the 2 by 2 Hamiltonian of Eq. 
Q can be interpreted as a pseudo-Zeeman interaction in val- 
ley space. While the spin part of the solution is simply |cr), 
the valley component of the eigenstates are the up and down 
eigenstates of the valley operator, ■ r. Then the four eigen- 
functions with energies -£'^,±,0- - Eo- ± AI-/2 are. 



1^, ±, cr) = ki, +) ® |cr> with ki, +) = 



ei't'KK' cos 



+e 



fKK' sm ■ 



(5) 

where the spinor in valley space for the |^, ±, cr) solution have 
Q'j_g. and (pj^j^, as inclination and azimuthal angles, respec- 
tively, in the valley Bloch sphere (see Fig|2]). The inclination 
angles of the (+) and (-) solutions with spin cr fulfill 
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ai^ = ai-ctan(A|,^,/5S) 



mod 71, 



(6a) 
(6b) 
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1^') l^'> 

FIG. 2. (a) Effective Zeeman-like field in the valley space Bloch 
sphere. The valley mixing is spin independent and it introduces the 
in-plane component; on the other hand the diamagnetic effect and 
the spin-orbit coupling enter in the out of plane component. The 
situation plotted corresponds to zero magnetic field, (b) Angles of 
the valley-space spinors corresponding to the two eigenstates with 
spin cr. 

which follow from Eq. ( |4bl l. 

We see from the expression Eqs. Q that when the spin- 
orbit coupling dominates, the z-components of the Zeeman- 
like valley field, 6^., have opposite sign for the two spin pro- 
jections. This situation is sketched in Fig|2|a), the unit vectors 
are drawn on a Bloch sphere in valley space for the B = 
case. Focusing on the two lowest (or the two highest) energy 
eigenfunctions, this implies that solutions for opposite spins 
have inverted weights in K and K' valleys. On the other hand, 
if l^orbl » |Afol, the sign of (5^ does not depend on cr. In the 
high parallel magnetic field limit |^| » a|.^, and the valley 
characteristic of the lowest energy solutions for different spin 
cr tends to be the same, i.e., they become valley polarized. 

It is worth noting that the Hamiltonians for opposite spin 
projections are related by the transformation: 

//^(B,) = Ti(4(-B,))*Ti, (7) 

which for B, - reduces to the time reversal symmetry re- 
lation. The last expression is useful for a direct construction 
of solutions for opposite spin projections (and opposite field 
Bj.) just by interchanging K and K' and conjugating the wave- 
functions. Explicitly, this means that for a known solution, 
l^icr, B-)}, there exists a related solution with the same en- 
ergy, |*F(cr, -B-)}, such that, 

mcr, B,)) = il^Ka \Kct) + iffK'a \K' Ct) , 

-BJ) = Ifrl'a \K^} + ^Kcr \K'^) > (8) 

which in terms of the inclination angles it reduces to recog- 

£ J: 

nizing that a\ a-iB.) - -{-B,). For B = the two degener- 
ated solutions 1^, d, T) and |^, d, i) fulfill Eq. To simplify 
the notation, we introduce in the following the doublet index, 
d = -H,-. The two states in the doublet d make a Kramers 
doublet because they are linked by time reversal symmetry. 



Each doublet can be regarded as a spin- ^ system, where un- 
like to spin-only quantum dots the the orbital parts of the two 
states are different. 

B. EfTective hyperflne interaction inside the doublets 

For brevity we omit in this section the dot index ^ = L, R, for 
instance the eigenstates d, cr) become |d, cr), and the below 
applies for both dots. The effective Hamiltonian seen by a 
confined electron in a quantum dot, as a result of the hyperfine 
interaction with the '^C nuclear spins, can be cast asj^ 

<F4i^'h''-' (9) 

,■={) 

where we have adopted an isotropic hyperfine field interaction 
(HFI), although there is some degree of anisotropy.'^The dy- 
namics of the electron spin is much faster than the precession 
time of the nuclear spins, T„,ic, and furthermore, each single- 
shot measurement in the RPE is realized over a time much 
shorter than r„„e. Therefore during each experiment the fields, 
h*'\ are constants given by the matrix elements of the hyper- 
fine interaction Hamiltonian, //hf = ^iso 2/ ' ~ R/)- 
[where Aiso is the hyperfine coupling constant and the summa- 
tion is taken over all the lattice sites in the quantum dot region, 
/, that have a '■'C atom with 1/ the spin- 1/2 vector operator of 
the nucleus located at R/] and tracing over the ensamble of 
nuclear spins. 

To obtain the RPE many realizations of the single-shot mea- 
surement are averaged over. The numbers /i*'*, with / = 0, 1, 2 
and i - X, y, z, follow Gaussian probability distributions with 
zero mean and the following variances, 

(10) 

where Nqd the number of atoms in the quantum dot and v is 
the abundance of '^C atoms in the dot. 

We now focus on dots that have both nonzero Akk' and 
nonzero Ajo splittings .t^^E^EHIl The estimation for variances 
in Eq. ( [TO] i are typically much smaller than the energy sep- 
aration, A = ^JiAxK')^ + (A.5o)^, between the d = -H and the 
d = - doublets. Due to this condition, the HFI mixing of elec- 
tron states belonging to different Kramers doublets can be ne- 
glected and therefore we only include matrix elements within 
the doublets d = +, 

(d,(r|H^ff |d,cr'). (11) 

This allows us to write an effective hyperfine interaction for 
the doublet d as: 

Htip^B'-a' + By,, (12) 

where the operators crjj, cr^, cr^, and are the Pauli matrices 
operating in the (d) Kramers doublet states and 
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^0 = {l^? cos ifiKK' - hf^ sin ipkr) (sin ffd.T - sin a^^y) , 

Bi = /.f COS (^ilZ^) + d cos - sin ^,,,) sin(^^^^l^) + d sin + /f > cos ^,,,) sin (^^1^5^), 
fi,^ = /.f cos (^^ilpi) + d (/.f ) COS - /f ) sin sin (^^ilpi) - d sin + cos sin (^^il^i), 



= /if + ^d cos - /if sin (sin aa.r + sin a^) . 



(13) 



These effective fields can be readily written as a function 
of the parameters of the Hamiltonian by using the expres- 
sions given in Eqs. (|6]). Furthermore, as discussed above, at 
zero magnetic field the inclination angles of the two solutions 
within the same doublet are related by time reversal symme- 
try, Q-djH-ad.i (they are opposite spinors in the valley Bloch 
sphere), and the expressions become simpler. 

The new fields also follow Gaussian probability distribu- 
tions with zero mean. The variances for the d = ± Kramers 
doublet's effective field (defined as cr^^ = ^(b^) ^, with 
a - 0,x, y, z) are, 



cr; 



cr 



= ^ (sin^ ad,T + sin^ ad,i - 2 sin Qfd,t sin ad,i) > (14a) 

(14b) 
).(14c) 



H 

2 2 

—- = — r- = 1 -H - smffd,! smad,i. 



d ^ 1/2 2 

= 1 H — (sin cKd.T + sin ad,i + 2 sin o-d.-f sin a^ i 



In Figj3]we plot the variances of the 3 components of the 
effective magnetic field as a function of the ratio Aso/^kk' for 
the case of zero magnetic field. The component proportional 
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FIG. 3. Effect of the valley mixing on the total root mean square 
of the effective hyperfine fields components (x, y and z directions) 
that act on the subspace of the two states in one of the &amers dou- 
blets. We also show the partial contributions that arises as a result of 
(gray solid line) the valley conserving components of the HFI, h*"', 
and (black dotted line) the non- valley conserving HFI components: 
the addition of the contributions generated by the h*'' and the h*^' 
components. 
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to the identity operator in the Kramers doublet only produces 
an energy shift and therefore it is irrelevant for the dynamics 
of the problem. The first important result is that the total vari- 
ance of the three relevant components are equal, moreover, 
they are the same in the two doublets. This means that the hy- 
perfine field within each doublet is statistically isotropic and 
therefore the dephasing of a given state in the doublet is sim- 
ilar to the dephasing of an electron spin confined in a GaAs 
quantum-dot. 

We also present the contribution of the valley-mixing terms 
(oc h*'' and oc h^^') and the valley-conserving terms (oc h*"*) of 
the effective HFI effective Hamiltonian of Eq.(|9]l. Both contri- 
butions are isotropic when the disorder dominates. This hap- 
pens because the eigenstates are valley spinors having equal 
weight in \K} and \K') (the spinors lie in the plane of the val- 
ley Bloch sphere) and therefore the two class of terms can mix 
those states. On the other hand, the states for Aso » Ai^k> are 

|-,T) = |/s:>®|T>, I- i> = |r)®|i), (15a) 

I+,T) = |^:'>®|t>, l+,i) = |^>®|i>. (15b) 

In this limit the valley conserving hyperfine operator, Toh'"* cr, 
contributes to B'^ since it provides diagonal matrix elements 
between |d, cr) and |d, cr') states; at the same time the operator 
is unable to flip the valley and it does not contribute to the 
effective fields and B^. On the other hand, the operators 
Tih*'* ■ cr and T2h*^' ■ cr both flip the valley index and therefore 
they do not contribute to B^ but they contribute to B^ and By. 



C. The (0, 2) and ( 1 , 1 ) eigenstates 

Here we give the two particle solutions for two electrons 
occupying the right dot and for one electron occupying each 
in absence of tunneling, i.e., in the high detuning limit valid 
for the preparation, evolution and measurement stages — see 
Fig[TJa). For simplicity, we number the four eigenstates and 
eigenenergies in the dot ^ = L, R as \^n) and E^, with n = 
1,2,3,4. We define a two particle Slater determinant built 
from states and ^'n' as 



(16) 



For two electrons in the right dot, the (0, 2) charge config- 
uration, only six independent states can be constructed using 
the four single-particle states |R«). These eigenstates and their 
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eigenenergies are, 
(0,2): 



(17) 



where n <n' and J/rr is the Coulomb repulsion energy. 

The eigenstates in the (1,1) charge configuration are the 
Slater determinants constructed from states |Rn) and |L«'). 
Therefore we have the sixteen states and energies: 



(1,1) 



feO , ^Ix' = fL + eR + £r„ + + f/LR, (18) 



with J/lr the Coulomb repulsion energy for electrons in dif- 
ferent dots. 

Energies in Eqs. ( fTS] ) and ( [T7] i are shifted by the gate- 
voltage controlled energies fL and fR. For simplicity we de- 
fine, 

e = - - t^RR + f/LR, (19a) 

£^1/ =\^^L + 3eR + f/LR + f/RR) , (19b) 

where the detuning, e, is the difference and Ej^y is the average 
between energies and Zi^''-^* (neglecting the part depend- 
ing on n and n'). Energies for (0, 2) and for (1, 1) now become: 



7(0,2) 



^«!n'^ ~ ^Av + e/2 + Er„ + Ei^„> 



e/2 + Er,, + Er 



(20a) 
(20b) 



In what follows the common energy shift, E^v, is omitted be- 
cause it is irrelevant for the dynamics of the RPE. 

In a clean system the single-particle eigenstates in both 
dots have the same spin and valley properties and therefore 
the left/right dot part of the (1,1) eigenstates can be separated 
from the «-space part (valley and spin spaces). This separabil- 
ity allows one to write the two-particle eigenstates as products 
of singlet and triplet functions in spin, valley and dot spaces, 
see for example the clean case in Ref. |44] Here, on the other 
hand, the valley disorder profile and therefore the single par- 
ticle eigenstates are different in the two dots and each two 
particle eigenstate — the Slater determinants given above — 
becomes an arbitrary linear combination of the dot/valley/spin 
tensor product states. This avoids identifying the eigenstates 
as singlet and triplet states in spin and valley spaces. However, 
we will identify singlet-like and triplet-like states according 
to how they behave when considering single-particle tunnel- 
ing mixing between (0,2) and (1, 1) states. In the next section 
we see that some particular linear combinations of these Slater 
determinants — within the subspaces generated by degenerated 
solutions — are the most physically relevant states. 



III. INTERDOT TUNNELING 



, i. These isolated-dot single-particle solutions become mixed 
when including tunneling between the dots. In order to pro- 
ceed, we introduce the identity and Pauli matrices in left/right 
dot space with / = 0, 1,2,3 (taking |L) and |R) as the -i-l 
and -1 eigenstates, respectively, of the operator's). Further, 
for referring specific blocks of the Hamiltonian we define, 
(T^ = 2-1 (cro + cr,), (t^ = l'' (cr,, - (r,), 'l = 2'' ('„ + fs), 
and 'r = 2-1(^0 -ft). 

We assume that the interdot tunneling preserves valley and 
spin degrees of freedom. This is because the tunneling am- 
plitude follows from the overlap between the quasi-bounded 
states in each dot. We take a tunneling energy, f, inde- 
pendent of the spin and valley that is being tunneling; this 
is a good approximation if the height of the confining bar- 
rier is much larger than all the energy scales in the DQD 
Hamiltonian.-^ Then the interdot tunneling Hamiltonian be- 
comes, Ht = -f^iCToTQ. The complete 8-dimensional single- 
particle Hamiltonian can now be compactly written as. 



^DQD - 'r^^O + ^L^^o ~ fftCToTO- 



(21) 



In a clean system the isolated-dot solutions of H^, |', d, cr), 
have the same spin and valley characteristics in the two dots 
and tunneling can therefore be described by four independent 
subsystems of 2 by 2 mixing |L, d, cr) with |R, d, cr), which al- 
lows one to build the two particle (0, 2) and (1,1) eigenstates 
as tensor products of triplet and singlet functions in left/right 
dot, valley, and spin spaces.l^SEHSESISI Moreover, the selec- 
tion rules for the mixing between (0, 2) states and (1,1) states 
follow correspondingly. There is only one avoided crossing 
produced by the interdot tunneling for each (0, 2) state with 
one (1, 1) state.l^ 

In the general disordered case considered here, the disorder 
in each dot can be different, and therefore the mixing becomes 
more complicated. We now proceed to the case of zero field 
(or B II z). According to Eq.([3]) the Hamiltonian within each 
dot does not mix opposite spin-z projections, and we have 

Hdqd = ^lo-^H\ + ^lo-iH]; + ^RO-fHf + ^Rcr^/Zf - f'lcroro. 

(22) 

Each valley operator have the valley solutions, |5^^,d^, 
given in Eq.(|5]l. Using as basis the eigenstates of each dot, 
we see that the interdot mixing matrix elements tT, = 
<L, dL, cr| Ht |R, dR, cr) are 



W,RdR ~ (^RdR,LdL) 



(23) 



COS COS 

2 



+1 sm cos 

2 



L,0- 



A. Disorder induced tunneling between states in different 
Kramers doublets 

The eigenfunctions of the quantum dot Hamiltonians, i/g 
(see Eq.([T])), are |ftd, cr) with ' = L,R, d = ±, and cr 



where <fRi - 1^*^/-<^^^/ ■ As in the clean case, the state |R, d, cr) 
can mix with the state |L, d, cr), but now also with the state 
|L, d, crj. This situation is sketched in Figj4|a). 

In order to derive properties of the tunneling amplitudes, 
we arrange the four quantities, 'l(]rjj//(-0, as a 2 by 2 matrix 
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(a) dot L 
^ ' E 



dotR 




(0,2) states 



FIG. 4. (a) Due to the disorder the valley mixing is different in the 
two dots, and each a state in Kramers doublet Ld is a different val- 
ley combination than the a state in doublet Rd. Nonzero tunneling 
amplitudes appears between a states in different Kramers doublets, 
(b) For a two-particle (0, 2) Slater determinant — which is an eigen- 
state in the high detuning limit — the interdot tunneling have nonzero 
matrix elements with four (1,1) Slater determinants. 



in doublet space, Dcrxs.- Using the unitary transformations, 
f/^, that diagonalize the valley operators i/^, one can readily 
show that D(r,LR is a unitary matrix. The same holds for the 
matrix D^-rl (= r)- This allow us to write. 



r -I- r 



LiM\ 



+ r 



(24) 



In addition, the orthogonality relation, a^^ + aL^ -n, com- 
bined (for zero magnetic field) with the time reversal symme- 
try relation, a* 



d,o- 



- lead to the following properties, 

d,(T 



'Ld.Rd 



Ld.Rd 



\f. \^-[f -V. 
\ Ld.Rd / \ Ld.Rd / 



(25) 



In what follows we use these expressions, and their relation 
to the parameters of the problem, to quantify the effect of the 
disorder in the probability experiment. 



B. Disorder induced avoided crossings in the two-particle 
spectrum 



In Eq.(|23) we have shown that the disorder induces single- 
particle interdot tunneling that does not conserve the Kramers 
doublet index. Before presenting the effect on the mixing be- 
tween the two particle states, we perform unitary transforma- 
tions in the doublet basis of each dot in order to work with a 
real tunneling matrix. This greatly simplifies the notation in 
the following. 

Instead of using as a basis the wavefunctions, |^, d, cr), pre- 
sented in Eq.Q, we multiply each of them by a phase factor. 



and the new matrix elements become, 

('Ld.Rd')„e„ = (fLd,Rd')old 



(27) 



There are infinite choices for making the new tunneling am- 
plitudes real numbers. We choose. 



/^^,. = Arg[(-rZ_,,_y, 



and then we have, 



fcdL = -^cos|, 



(28a) 
(28b) 
(28c) 
(28d) 

(29a) 
(29b) 



with ir\(j the angle between the vectors in valley space: ^ 
and The amplitudes (^RdL^') follow straightforwardly 
by conjugating the matrix elements above. 

A physical intuitive (and equivalent) approach for writing 
the tunneling amplitudes as real numbers follows from: (i) 
perform a unitary transformation in both dots making the z- 
axis coincide with one of the valley vectors, for example g^; 
(ii) rotate both dots around the new z-axis in order to take the 
other valley vector, i.e., 5^, to the xz-plane where its eigen- 
vectors can be written as real valley spinors fully described 
by the angle, rja-, between the vectors gir; (iii) leave the valley 
basis fixed in the first dot (^ = R) and perform a unitary trans- 
formation in the other dot (^ = L) making the real eigenvectors 
the new basis; (iv) the resulting tunneling amplitudes become 
real. 

From hereon we use the new single particle eigenstates of 
the isolated dots presented in Eq.(|26| for the derivations in- 
volving two particle states. 



C. Landau-Zener physics for different kind of crossings 



We now apply the tunneling Hamiltonian, its two particle 
version [H^^ = li ® {H^^)2 + {h}^)i ® I2], to the (0, 2) eigen- 
states. 



|R,d,o- \ 


y \r 1 


|R,d,,o-i/ - , 


/ i ['Ld'.Rd] 1 


d 


'=± 



R,d,o- \ 
L.d'.o-]/ 



+ f 



Ld'.Rd 



iL.d'.o- \ 
R,d,,<r|/ 



(30) 



|^,d,cr)„,„ = e''^-.^ I^,d,(r)„id, 



(26) 



This result is sketched in Fig|4|b) where we show that, by 
virtue of the disorder, each (0,2) Slater determinant have 
nonzero matrix elements with four , instead of two, (1,1) 
Slater determinants. 

We now focus on the disorder induced Landau-Zener pro- 
cesses that arise due to the additional mixing terms. For sim- 
plicity, in this first approximation to the problem, we assume 
that the energies, £/, of the four (1,1) Slater determinants (we 
label them / = 1, 2, 3, 4) are well separated on the scale of the 
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FIG. 5. (a) Simplest case for an avoided crossing as a result of the mixing between a (0, 2) state (|a>) and an (1, 1) state (\b)). A detuning pulse 
take the system, initialized in \a), from positive to negative detuning. The Landau-Zener probability to follow the avoided crossing is, P^^'. (b) 
Probability, Pa (in a DCAC) as a function of /'jj (in a disorder induced avoided crossing, i.e., a DFAC), for different ratios between the gaps, 
P = (Ajj/Ajd)^, assuming that both processes are realized at the same detuning speed. One sees that /?, that quantifies the disorder, must be very 
small for the probabilities Pjg to be neglected and for simultaneously having that Pa is close to 1 (adiabatic process), (c) Color plot of /? as a 
function of the parameters of the problem for B = 0. We see that moderated values of the valley mixing parameters can generate non-negligible 
values of p and therefore the Landau-Zener physics in the disorder induced avoided crossings must be in general taken into account. 



tunneling t and therefore we can write four separated effective 
Hamiltonians describing the mixing with the (0, 2) state as. 



H, = 



2 

-tl 



-tl 

+ Ei 



(31) 



At the crossing energies the mixing due to interdot tunnel- 
ing is maximum and avoided crossings with gaps A/ = 2|f/| 
appear. By inspecting Eqs.(|29]l and ( |30l ) we readily identify 
two different avoided crossings types; (a) the doublet conserv- 
ing avoided crossings (DCACs), with a gap value of, AJj"^ = 
2|f cos(77o-/2)|; and (b) the doublet-flipping avoided crossings 
(DFACs), with a gap value of, A|[- = 2|f sin(77<r/2)|. The latter 
type, DFACs, is the one induced by disorder because when the 
two valley vectors are the same, as in a clean case, we have 
77o- = and A^j^O. 

As described in the Introduction and in Figjl] detuning 
changes are intrinsic to the RPE single-shot cycle: the sep- 
aration and joining stages. The gap is an important parameter 
for determining the probability of realizing an state conversion 
when a change of detuning is applied and the system passes 
through an avoided crossing. In a given stage of the problem, 
h (with h - s,j, separation and joining stages, respectively), 
the Landau-Zener probability of realizing the state conversion 
in an avoided crossing, I, is. 



P*''' = exp 



IttAj 



(32a) 
(32b) 



The probability P^''^ grows the slower it is the rate of change 
of detuning at the crossing /, v^'"', and the bigger is the gap 



of the avoided crossing. A;. For convenience, we have de- 
fined * as the probability to remain in the original state, 
this probability approaches zero the more adiabatically the 
Landau-Zener process is realized. 

Because the gaps in the exponents of the Landau-Zener for- 
mula appears to the second power we define, 



f^^(^) =tan^(W2), 



(33) 



This number allows us to quantify the effect of the valley mix- 
ing for the Landau-Zener physics involved in the experiment. 
How sensitive is y6°" with the parameters in the double dot is an 
important question. Since the angle between the valley fields 
is given by. 



Tja- — arccos 



AIM 



(34) 



by virtue of the half angle relations, we find that /3'^ is given 
by the quotient. 



ALa;^-AL^,A^^, cos <p^i^-6];6l 
ALAR + 



(35) 



When B = 0, there is no dependence with cr in the angle 
T/o- and consequently, nor in the tunnelings amplitudes and 
the gaps. In what follows, we assume that the lengths of 
the two dots are the same so that the spin-orbit splitting for 
the lowest lying bound states is the same in both dots, and 
(5)^ = 5)^ - crAso- However, it must be noted that if the dots 
have different lengths leading to different spin-orbit splittings, 
then the effect of disorder in the Landau-Zener physics can 
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be even more drastic than what we show below, because the 
angle rja- becomes more sensitive to the in-plane valley fields 
diff'erences. Our formulas in Eq.([35| and ([34]) are general and 
can be used for such cases. 

Figure |5j a) shows a sketch of an avoided crossing between 
a (0,2) state and a (1,1) state reflecting the Hamiltonian pre- 
sented in Eq.([3T|l. In the example, the avoided crossing mixes 
an (0, 2) state, |a>, and with an state (1, 1), \b}. Note that the \b) 
Slater determinant differs from the |a) one in that a cr single 
particle right-dot eigenstate with double index d is replaced 
by a cr left-dot eigenstate with the doublet index d'; thus the 
mixing is due to the f^^, tunneling amplitude. The system is 

prepared, at positive detuning, in the state |*Fp^ = \a) and then 
the detuning is swept to negative values. The Landau-Zener 
process governs the probabilities that the final measured state, 
\^f), reflects a state conversion, \^'f) cc \b) (with Pdd' as given 
in Eq.(|32]i), or it remains in the original state, |^f ) oc \a) (with 

^dd' = 1 - ^'dd')- 

Figure |5jb) shows the Landau-Zener probabilities (see 
Eq.([32]i) in a DCAC versus the probability in a DFAC linked 
by different values of /3 assuming both Landau-Zener pro- 
cesses were performed at the same detuning rate. The result 
shows that only for /3 <K 0.01 one could ignore the disorder 
induced avoided crossings assuming that P^j x and at the 
same time assume that the doublet conserving crossings pro- 
vide adiabatic state conversions, i.e., Pm ~ 1. In panel (c) we 
present a map of /Jo- (- /Sg- - /3), from Eq.((35]l, for diff'erent 
phase differences between the valley mixing terms and as a 
function of their strengths, A|,^,, normalized to the spin-orbit 
coupling splitting. The plot allows to see that p can be above 
0.01 — i.e., the Landau-Zener physics in the DFACs cannot be 
neglected — in a broad region in parameter space, even in sit- 
uations in which the spin-orbit coupling dominates in both 
quantum dots. 



D. Full spectrum, classification of the avoided crossings 

As shown above for the zero field case, the Landau-Zener 
physics in doublet-flipping avoided crossings can be impor- 
tant for the return probability experiment and therefore their 
effects must be studied. If the tunneling energy, t, is larger 
than the splitting energies in the dots, A^, the mixing of each 
(0,2) states with the four (1,1) states given in Eq.([30| must 
be consider simultaneously. No simple Landau-Zener physics 
concepts can be applied because the mixing cannot be decom- 
posed into two-level avoided crossings. In what follows, we 
study the problem assuming that t is smaller than the A^ ener- 
gies and that the valley mixing is sufficiently different in the 
two quantum dots so that |A^ - A'^l > t. In this limit, the 
22 levels two-particle spectrum (describing the mixing of the 
six (0,2) states with the sixteen (1, 1) states) can be decom- 
posed into two-level avoided crossings. The LZ formula can 
be used to describe all the processes in the different avoided 
crossings happening as the separation and the joining stages 
of the PRE are performed. The qualitative conclusions drawn 
for the low t picture can be extended to the general t case. 




FIG. 6. Avoided crossings classification in the two particle spectrum 
for zero magnetic field reflecting the mixing between the six (0, 2) 
states and the sixteen (1,1) states as a function of detuning. The 
avoided crossings III, III', II and IV are induced by disorder, they rely 
on tunneling matrix elements that do not conserve doublet index. 

We also note that for zero magnetic field the tunneling ampli- 
tudes given in Eq.(|29|l can be chosen to be independent of cr, 
therefore in what follows we omit the superindex cr and set 

fLd'.Rd = fLd'.Rd- 

The small t assumption allows us to introduce a classifica- 
tion of the avoided crossings. In Fig|6]we show such a classi- 
fication in the two particle spectrum as a function of detuning, 
obtained by taking the two particle version of the Hamiltonian 
of Eq.(|2T]i and using the Slater determinants constructed from 
the valley and spin single particle states, l^rcr) - 1^) ® |t) ® |cr), 
i.e., the six independent (0,2) states with (T'o-')i^(Tcr), 

plus the sixteen (1, 1) Slater determinants, 1^^^^^,^ The avoided 
crossings presented in the figure are studied in detail below, 
their associated detuning and energy values are given in Table 

m 

From Figure [6] it is clear that in the high detuning limit the 
sixteen (1,1) states are grouped into four subsets of four states 
given by which doublets are occupied in the left and the right 
dot, dL and dR, respectively. We denote those four subsets as, 

(LdL,RdR), with states, 1^^^^/). and energies, 

^dJU "'^/^ + dLAV2 + dRA'^/2. (36) 

Similarly, the energies of the (0, 2) states in the high detuning 
limit are, 

<d'k =-«/2 + diA^/2 + d2A'^/2. (37) 

For di = -d2 the energy level, -e/2, is fourfold degenerated 
since it is associated with the four states, In^"^, V For future 
reference we denote this subset as (R+, R-). 

We first focus on the avoided crossings involving the two 
(0, 2) non-degenerated states (with high detuning energies in 
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Eq.((37]i taking di = d2 = +), namely, the ground state and the 
highest excited (0,2) states. 



Rdi/ ■ 



Because the two electrons occupy the same Kramers doublet 
d the total spin is zero. We use the letter S for these singlet- 
like solutions even though these states are not singlets in spin 
space: the spin part of the two particle state cann ot be sep- 
arated even in the absence of valley mixingl^SEl Following 
Eq.([30ll and applying the low t picture, we can treat separately 
the mixing due to the tunneling Hamiltonian between the 



states IS'I,"'^'^ and (1,1) states in different levels (LdL,RdR). 



Since the single particle tunneling is unable to change two 
single particle states simultaneously, the state jS'ly'^*^ does not 
mix with the states in subsets (LdL,Rd). 

Each doublet conserving tunneling amplitude, fLd,Rdj mixes 
the state with the (1,1) combination. 



.(1,1) 



\ = J_/'|RdT\ 



LdT 
Rdi 



(39) 



of the subset (Ld, Rd). This defines the I (d = -) and F (d - 
+) doublet-conserving avoided crossings shown in Figure |6] 
In this case the matrix element of the tunneling operator is 
V2fLd,Rd and therefore the gaps are Ai = Ap = 2V2I cos = . 

Similarly, each tunneling amplitude, f^a Rd' mixes the state 
I'^dd^') ^^'■^ '■^^ combination. 



.(U)\ ^ 

dd / 



1 /iRdT 

V2 V"^ 



LdT 
Rdi 



(40) 



of the subset (Ld, Rd). This mixing generates the two doublet- 
flipping avoided crossings. 111 (d - -) and 111' (d = +), with 



2V2 



sm : 



gaps Am = ^uy - - 21 

Apart from the factor V2 in the gaps, the avoided crossings 
I, F, 111 and 111' share another important feature. Because in 
each crossing the single (0,2) state, j'Sly'^*), crosses four de- 
generated (1,1) states (and mixes with only one combination 
of them) three other (1, 1) states of the (Ld',Rd) subset are 



TABLE I. Detuning and energy values for the avoided crossings 
presented in Fig|6] The detuning position of the avoided crossings 
are obtained by equating the high detuning energies of Eq.l|36[) and 
Eq.l[37j. The Kramers doublet index conservation (or not) of the 
tunneling that generates the mixing is also stated. 

Ci,C2 £ci =gC2 Ec^ Ec2 Type 



III, ir - 7 (a^ + h^) 

IIF, II i ^Ai' + A'^) 

I, IV' \ - h^) 

r, IV \ [h^ - A'-) 

V 



I (a^ - 2,h^) \ (a^- + A"*) DFAC 

\ (3A'* - A'-) -\ + h^) DFAC 

- i (a^- + \ [h^ - A^-) DCAC 

\ (a^ + i (a^- - A'*) DCAC 

- DCAC 



(a) (GA M = (!'-)■ (i',++),(ni,+-) or (iir, -+) 



(38) 2V2|?,,^^,jJ 




(C) (i) A" > A^ (ii) A" = A^ (iii) A'' < A^ 









7>i 







I R+o-y I L+o-/ ^ I R+(T' 

FIG. 7. States involved in the avoided crossings presented in Fig|6] 
(a) Avoided crossings G = 1,1', III or III', the last two are induced 
by disorder. As a single (0,2) state crosses four (1, 1) states there 
are three (1,1) blocked states — not mixed by interdot tunneling — 
in close analogy with spin-only double dots, (b) No blocked states 
appear in the disorder-induced avoided crossings II and IF, four (0, 2) 
states cross four (1,1) states, (c) The doublet conserving avoided 
crossings IV and IV' also do not present blocked states. According 
to the values A^ and A"* the detuning dependence changes which 
lead to different outcomes in the RPE; we refer to the grey arrows 
in Sec jlVB 2| when showing the two possible situations available if 
|5f;"'^ is prepared. When A'^ = A'' (extremely unlike in a disordered 
system) the (L+,R-) and (L-,R+) subsets become degenerated and 
IV and IV' would collapse in a single crossing, V. 



unafi^ected by the interdot tunneling. These states are. 



y +i,d'd/ - 
U(l,l) \ ^ 

r -l.d'd/ - 



d'T\ 
A} I-- 

d'i\ 
All-' 



Ld'1 
RdT , 

Ld'J 
Rdi , 

1 



|t(i.i)\ = _L /iRdT \ 

l-'o.d'd/ - ^[^ l|Ld'i/ 



V2VI' 



Ld'T 
Rdi 



(41a) 
(41b) 
(41c) 



with energies growing linear with detuning following 
Eq.([36|. We denote them with T in direct analogy with the 
spin-triplet states of the (1,1) configuration in spin-only dou- 
ble dots; in such a situation the interdot tunneling only af- 
fects the spin-singlets (0,2) and (1,1) leaving the spin triplet 
states unaffect ed. In the language of Pauli-blockade these are 
blocked states. 02211! We will refer to them as triplet-like states 
despite the fact that |7"q'(J!^) is not a spin-triplet. 

We now focus on the crossings involving the remaining 
(0,2) states, i.e., states of the (R+,R-) subset. Again, the 
small t picture allows us to study separately the four cross- 
ings with the (LdL,RdR) subsets of the (1, 1) configuration. 
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By inspecting the action of the interdot tunneling operator on 
an (0, 2) eigenstate given in Eq.(|29| one can show that all the 
crossings become avoided crossings and no blocked states re- 
main. This is easily seen by noting that any (1,1) state 
is mixed in an effective 2 by 2 system- 



gaps. 



neling amplitude, t,j dj — with the (0,2) state C^"'^,). 



-due to interdot tun- 

,Rd^— wiu. u.c vu,z.; SLaLC The 

singlet-like and triplet-like basis defined in Eq.([39]l, (|40]l, ( |4T] i 
facilitates the calculation of the hyperfine dynamics. For the 
avoided crossings here, we have 



\ o-,,dLdR| T I o-f.dRdR 
\ dLdRl T I dRdR/ 



LdL,RdR' 



(42a) 
(42b) 



Ai = Ai- =2V2|fL-,R- 
Ani = Am' = ^J/3^l, 

Aiv = Aiv = -^^i' 



A„, 




(44a) 
(44b) 

(44c) 
(44d) 



In the next section we make use of these relations for calcu- 
lating the outcomes of the multiple Landau-Zener processes 
in the separation and joining stages. 



IV. RETURN PROBABILITY, METHOD AND RESULTS 



with cTf - - 1 , 0, 1 and, 

1.^(0,2) \ ^ |RdT\ U(0,2) \ ^ |Rdi\ 

r +l,dd/ - iRdT/ ' r-l,dd/ - iRdi/' '^^^> 

|7^(0,2)\ ^ J_ /|RdT\ _ |RdT\\ L(0,2)\ ^ J_ (\Rdt\ |RdT\\ 

\ om/- ^Wm/ I'^dd /- ^Wm/ i^'^i/)' 

i.e., we have defined triplet-like and singlet-like states inside 
the fourfold degenerated subset (R+, R-). 

In particular, the doublet conserving tunneling amplitudes 
generate avoiding crossings with (1, 1) states in the (L+,R-) 
and (L-,R-i-) subsets; in Fig|6]we have label these avoided 
crossings as IV and IV', respectively. Their gaps are given 
by Aiv = Ajy = 2|fcos jj. Conversely, the doublet-flipping 
avoided crossings II and IF have gaps. An - Aiy - 2 \tsm j\, 
and are associated with the crossings between (0, 2) states be- 
longing to the (R-i-,R-) subset and (1, 1) states belonging to 
the subsets (L-, R-) and (L-H, R+), respectively. 

In Fig|7]we present an detailed sketch of all the mentioned 
crossings that include the relevant states. Panel (c) reflects the 
fact, as it is shown in Table |l] that the relative positions of 
crossings IV and IV' depends on the difference between A** 
and A^. This happens because the avoided crossing involving 
the (0,2) subset (R-I-, R-) is IV when crossing (L-I-, R-) and 
IV' when crossing (L-, R-I-). As it is shown below this leads 
to different outcomes of the return probability experiment if 
excited (0,2) states are prepared. Finally, in the central panel 
of Fig|7|c) we show that IV and IV' collapse into a single 
avoided crossing, V, when A** = A^. We do not consider this 
avoided crossing further because such an special case is not 
representative for a disordered system. 

It is worth to note that in the low t picture the knowledge 
of y6 allows one to quantify all gaps in the 22 by 22 system as 
a function of a single one. We choose as reference the valley 
conserving avoided crossing I. In a clean system the energy 
gap, 2 V2 is common to the six existing avoided crossings, 
namely, I, I' and the four crossings associated with case V"^"* 
Here we have eight avoided crossings and four distinct energy 



A. Obtaining the return probability 

Having identified the complexity added by the presence of 
valley mixing in the double dot we now develop a scheme for 
computing the outcome of the return probability experiment, 
including all stages. The path we take is designed for gaining 
knowledge about the physics of the experiment with a degree 
of generality, any quantitative study designed for fitting a par- 
ticular experiment must include further details as the actual 
pulse shapes and the dependence of the parameters with de- 
tuning. A full numeric time-dependent approach would be 
necessary if the avoided crossings overlap each other avoid- 
ing the application of the Landau-Zener formula separately 
for each crossing. 

The most important information that the experiment pro- 
vides is the dephasing time, Tj. It is also important the shape 
of P{ts ) and in particular its saturation value, 

Po. = Hm P(t,). (45) 

Our theoretical study aims for a qualitative understanding of 
which saturation values can be expected due to the hyperfine 
field induced dephasing in a disordered nanotube double dot. 
We show that the avoided crossing induced by disorder can 
also affect the return probability value even in the absence of 
dephasing. If the probability, 

Po = hm P(t,), (46) 

is smaller than 1 it means that the application of the separation 
and joining stages — without waiting at the evolution stage — is 
reducing the return probability due to the additional Landau- 
Zener processes. As it is discussed below, this statement is 
relevant for prepared states different than the ground state. 
Figure [8] sketches two situations studied below, preparing the 
ground state, panel (a), and preparing an excited state, panel 
(b). 

We focus now on the preparation stage. Any given 
(0,2) prepared state can be decomposed into the six (0,2) 
eigenstates of Eq.([T7]i, that we label |Z, (0,2)), as, ITp) = 
YXi (0' 2)). To obtain the return probability one computes 
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FIG. 8. Single-shot return probability experiment: separation ("s"), evolution ("e", during a time T5) and joining ("j") stages. Avoiding 
crossings (A) are doublet-conserving and, in the situation we study, their gaps are much bigger than the others; the state conversions in 
crossings (A) are almost adiabatic. (a) The prepared state, Tp, is the (0,2) ground state and the separation stage involves the avoided crossing 
I which is type (A); at the beginning of the evolution stage the system is in a (1,1) state. Due to the hyperfine dephasing during the evolution 
stage (e, Tj ), states blocked for crossing I are populated. Those blocked states can still return to (0, 2) at the joining stage due to the disorder- 
induced avoiding crossing II. The return probability thus depends of the detuning value at the measuring stage, see that P' t f", see Sec IV B 
(b) The prepared state is the (0,2) highest excited state. At the separation stage the system passes through multiple Landau-Zener process, 
moreover, the first one is the disorder-induced avoided crossing, III'. The system is not guaranteed to be in a (1,1) state at the beginning of the 
evolution stage, we study in detail this case in Sec |IVB 2] 



the average, over an ensemble of random hyperfine fields, of 
the single-shot measurement's outcomes; the final result does 
not depend on the phases of coefficients 

P^„(r,)^Yj\a,fP'(^-^^^ (47) 

where P;(ts) is the return probability found if the (0,2) eigen- 
state |/, (0, 2)) is prepared. We then study separately the RPE 
assuming that a given eigenstate is prepared. Equation ( |47j i 
can also be applied to prepared mixed states, e.g. a ther- 
mal mixture, by replacing with the probability weight of 
the state \l, (0, 2)) in the prepared mixture. However here we 
will assume that the < A**^^ (for typical systems and 
a'" are a few hundreds of yueV and the temperature is below 
lOOmK, (k^T < 9jueV)) and therefore individual eigenstates 
can be prepared.'^ 

When preparing (0,2) excited states, see Fig|8jb), the sys- 
tem must pass through multiple Landau-Zener processes at the 
separation stage and similarly at the joining stage. As shown 
in Fig|8ja), this is not the case if the ground state is prepared. 
In order to solve this problem for a general prepared state, we 
profit from the low t picture and from the fact that the ex- 
periment is performed over an ensemble of hyperfine realiza- 
tions. This allows us to obtain the final result by analyzing the 
Landau-Zener processes from a probabilistic perspective (free 
of interference effects) and combining those (both separation 
and joining stages) with the averaged behavior of the dephas- 
ing due to hyperfine at the evolution stage. It is important 
to keep in mind that t is still much bigger than the hyperfine 
characteristic energy cth of Eq.( fTO|), the duration of the sep- 
aration and joining stages (see Figmc)) is shorter than ti/crii 
and therefore it is still a good approximation to consider the 
effect of the hyperfine interaction only at the evolution stage. 

When applying the probabilistic approach, we are assum- 



ing that the quantum interference effects, that arise only in 
Landau-Zener loops (when the system can arrive to a given 
state through different paths), average out. In order to explain 
the approximation we focus on the example shown in Fig|8jb) 
in which the highest excited (0,2) state is prepared. At the sep- 
aration stage, the prepared state first visit the avoided crossing 
III' and then the crossing 11' (note that I' and IV' are doublet 
conserving avoided crossings that in the figure are assumed as 
adiabatic processes). The energy difference between the two 
branches of the loop is given by A''. If, during the separa- 
tion detuning pulse, the time spent inside the loop, Tj, varies 
ATi, the phase difference that controls the interference after 
the crossing IF changes in A0 x A^ATi/H. The interference 
averages out if the phase gets randomized in a range A<p > 2n. 

Taking the typical experimental values for A^, the random- 
ization happens if ATi varies in a range of at least 15ps. Such 
ATi variations are expected to be present in the experiment 
when combining the errors in preparing the system at point 
"p" (see Fig{T^) with the imperfections of the detuning rising 
pulse from "p" to "e". If, on the other hand, ATi < lOps, the 
phase changes are not big enough and the interference effects 
would persist.'^ In such cases the single-shot cycles can be 
designed ad hoc with different times T/ in order to get rid of 
the interference in the Landau-Zener loop and the treatment 
we do here would still be valid.^The main conclusions of the 
paper regarding the preparation of excited states is that the ex- 
periment no longer responds to the standard return probability 
experiment, leading to Pq < 1 ■ Those conclusions apply to the 
general case: we have performed full time dependent simula- 
tions (that contain the interference effects) and, if they are not 
tuned to sit in a destructive interference condition, they share 
this feature. 

By working in the limit that the quantum interference in the 
Landau-Zener loops can be neglected, we manage to analyze 
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in a simpler frame the general features that arise when prepar- 
ing excited states. We define the matrices M**' for the sepa- 
ration stage, M^''\ts) for the evolution stage, and M® for the 
joining stages. We make these matrices explicit below. They 
operate on the space of the 22 eigenstates with zero-tunneling, 
i.e., the six (0,2) states (that can be prepared) and the sixteen 
(1,1) states. The result of operating with these matrices on a 
given vector gives the probabilities of finding the system in 
a new state after performing the stage. Therefore the vectors 
and matrices here are real positive numbers; in what follows, 
for each regular quantum state we associate a real probability 
vector using the modified ket and bra notation. 



(48) 



We have introduced the ~ symbol to clearly distinguish the 
probabilistic picture from usual quantum mechanics. This 
non-quantum mechanical bra and ket notation greatly simpli- 
fies the writing of the matrix elements of the probability ma- 
trices and their operation on probability vectors. 

Following Eq.pT) |fl;p are the nonzero components of |*l'p), 
the probability vector associated with the prepared state; the 
components lie in the (0,2) block of the vector Since in what 
follows we study the RPE after preparation of specific (0,2) 

eigenstates then a single component of j^Fp^ is nonzero and 
equal to 1. We apply the separation stage and obtain the real 
vector. 



(49) 



that contains the probabilities of finding the system in each of 
the 22 states at the beginning of the evolution stage. The ma- 
trix M^'^Xts) is applied next, it contains the transitions prob- 
abilities between the different states due to the action of the 
hyperfine field during a time r,. The average over the hy- 
perfine field ensemble is already contained in M^^\ts). The 
probabilities after the evolution stage are. 



K(T.))=M('='(T,)k^ 



(50) 



And finally we apply the joining stage matrix to obtain the 
final probability vector. 



(51) 



The probabihty to measure the system in an (0,2) charge state 
is obtained simply by summing over the six components asso- 
ciated with (0,2) states. We define the (0,2) real vector proyec- 
tor, |(0, 2)^ having the six (0,2) components equal to 1 and 
zero for all the other ones. The return probability then be- 
comes, 

P{T,)^{(Oa)\¥r,). (52) 



1. Separation and joining stages 

The states involved in each avoided crossing follow from 
the analysis presented in Sec 111 D From Table|l]and Fig|6]we 



recognize four different detuning values, s\, s\i, em and siv 
in which two avoided crossings operate. As an example, at 
e\ the action of the avoided crossings 1 and IV are comprised 
in the matrix, , where (s) refers to the separation stage. 
Since the crossing I is a Landau-Zener process between states 



1^(0,2)^ and jS'':'!'*^ we have. 



,(0,2) 



,(0,2) 



,(0,2) 



,(1,1) 



^(1,1) 



,(1,1) 



:(i,i: 



,(0,2) 



Pf, (53a) 



□(s) 



(53b) 



Similarly, given the mixing at crossing IV described by 
Eq.(|42]i, we have. 



,(0,2) 

n(03T 
o-,,-+ 

(iTT) 

-+ 



,(0,2) 



^(U) 



,(1,1) 



-(s) 
IV' 



M 



(s) 



n(l,l) 



n(0,2) 



,(0,2) 



n(0,2) 



n(l,l) 



•,(1,1) 

- o-,,-4 



IV' 



^(0,2) 



n(0,2) 



M 



(s) 



n(l,l) 



5(s) 
IV ' 



(54a) 
(54b) 
(54c) 
(54d) 

with cr, = -1,0, 1. All the remaining matrix elements of Mg^^ 
are trivial, being the non-diagonal ones and 1 the diago- 
nals ones. We construct in the same way the remaining three 
matrices of the separation stage, M^^, Mg^'p and M^^J. The 
matrices for the joining stage follow from replacing the LZ- 
probabilities with the ones corresponding to the joining stage 
P®. These probabilities can be different because, even though 
the gaps of the avoiding crossings are the same, we allow for 
a difference in the detuning rate of change between the sepa- 
ration and the joining stages. We define the ratio between the 
speeds as. 



v(i) 



(55) 



For simplicity we assume that the detuning speeds v^^' and 
v'*J' do not depend on detuning fast enough to be considered 
different at the four relevant detunings (i.e., in the detuning 
range [em, en], which from Table |l] is [-(A^ + ^^)I2,{A^ + 
A^)/2]). One can write all the Landau-Zener probabilities as a 
function of a single probability, for instance P^^\ once k and yS 
are known. We do so by combining the Landau-Zener formula 
of Eq.(32|, the relations fixed by y6 between all the gaps in 
Eq.(|^ I and the definition of k given in Eq.(|55]). 

The full separation (joining) stage matrix is then given by 
the successive application of the matrices above in decreasing 
(increasing) detuning order. 



(s) 



{mIm®m®m^I 



for > ^^ 
for hy <A^' 

for A^ > A"^ 
for A^ < A"^ 



(56a) 
(56b) 



14 



2. Evolution stage, hyperfine interaction 



As discussed in Sec |IIB[ during each single-shot cycle the 
hyperfine fields can be considered fixed. Because » cth, 
f - L,R, the evolution of an electron is assumed to be re- 
stricted to the Kramers doublet it occupies, i.e., the dynam- 
ics follows Hamiltonians or //^p" as the one given in 
Eq.([T2|i. The phase prefactors incorporated in the unitary 
transformation in Eq.p6l) modify the effective fields presented 
in Eq.([T3]l. This does not introduce any complication because 
the transformed Hamiltonians have the form of Eq.([T2]l with 
transformed hyperfine field components that follow the same 
zero mean Gaussian distributions of Eq.([T4]); i.e., the effective 
hyperfine field seen in any given Kramers doublet retains its 
statistical isotropic property. 

Therefore, the time evolution operator (taken from fo - 0) 
for an electron in the doublet (^d^) is. 



Mf^^(f) = exp(-ii/^'jff) 



cos(w^''f f)crf ' - i sinCw^^-'f) cr' 



(57) 



where cr^'^f s (cr^''^, crf^'' , crf^) are the Pauli matrices operating 
on the pseudo-spin of the (^d^) Kramers doublet, B^'^^n^'^^ = 
B^'^f is the effective Hyperfine field vector with n^'^f being the 
unit vector that defines the direction of the field, and w^''^ - 
B^'^'^in, the precession frequency. 

In order to compute the dynamics for any (1,1) or (0,2) state 
we use the Slater determinants' basis constructed with eigen- 
states of the two quantum dots, j^dfcr^^. The time evolution 
of each Slater determinant is dictated by. 



U{t) 



{t)xui%At) 



|f'df-o-2 



u 



fdf;f'df, 



(0 



If'df'a"! 
fdfO-i 



(58) 



From Eq.(|50| the matrix element Pijir,) - ^<E)/| M'-''\t,,) jO/^ 
is the probability (averaged over all the single-shot cycles) of 
finding the system after a time Tj in the state \(^f^ given that 
the state at the beginning of the evolution stage is From 
the doublet conservation dynamics dictated by Eq.(|57]l, it fol- 
lows that nonzero matrix elements arise only if the occupied 
Kramers doublets — for example ^idi and ^2'^! — in 1*1*/) and in 
\^f^ are the same. 

The probabilities Pij(Ts) are obtained by applying the time 
evolution operator of Eq.(|58|l to the two-particle state |<1),), and 

averaging [/(r?) over the hyperfine field realiza- 
tions that are independent in different dots and doublets. As 



shown in Sec 11 B for zero magnetic field the effective hyper- 



fine field components in a given dot follow zero mean Gaus- 
sian distributions with the same variances irrespective of the 
Kramers doublet, i.e., cr^ s cr^df,; with j - x,y,z and d^ - ±. 
All the probabilities can be expressed in terms of only two 



ensamble averages per dot,l2ll22l 



1 + 



1 -4 



hyperfine 



exp 



-2 



hyperfine 



1 -4 



" ^' s 



exp 



(59) 



(60) 



Note that the average over the hyperfine fields of 
(n^'^^) sin^((yf''''T.5) — where n^^'' is a cartesian component of 
the unit vector n^''^ — is independent of j because of the 
isotropy, in average, of the effective hyperfine field. Due to 
the zero mean of the Gaussian distributions that govern all 
components of B^''^ , terms having odd powers in any given 
component n^^'^ average out being absent in Pij{Ts). 

The nonzero probability elements for the evolution depends 
on T, as sums of products of the above functions. Starting with 
spin polarized states the matrix elements are. 



±l,did2 
- Tl,d|d2 
0,d|d2 
did2 



(e) 



±l,did2 
±l,did2 

±l,did2 

WZNr) 
±l,did2 



^ ifc + fl')(fc + ff\ (61a) 

= 4/|'/|% (61b) 

-2fl'fl^-+fl'f^^+f^'fl\ (61c) 

-2fl^fl'-+fl'fc+fcfl'-^ (61d) 



while if the evolution starts in singlet like states we have. 



-(Nl,Nr) 
did2 

0,did2 

r(ivDv«) 
±,did2 



ANl,Nr) 
did2 

'd,d2 

did. 



= ./^'4^+3/|'/f, (62a) 

-fcfi+fsft-+^fsff^ (62b) 
^fl'&+fl'&+2fp&, (62c) 



where for (Nl,Nr) = (1,1) there are four possible dida cases 
with ^1 = L and ^2 = R, while for (A^^, A^^) - (0, 2) only the case 
did2 = -I — is valid having ^1 =^2 = R- The matrix elements of 
M^^\ts) for the case that the evolution starts in a To-like func- 
tion can be obtained from Eq.(|62]) by making the replacement 
To <-> 5 . Two nonzero matrix elements remain to be defined, 
those for (0,2) singlet-like states in the same Kramers doublet. 
Since these are nondegenerated states that are well separated 
in energy from the remaining 21 states, we have. 



,(0,2) 
'dd 



(e) 



(0,2) 
dd 



= 1. 



(63) 



B. Results, zero magnetic field 

The separation and joining stages produce different out- 
comes for singlet-like and triplet-like states even if they be- 
long to the same subset (^idi,^2d2). This can be seen in 
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Fig|7|a) for crossings I, F, III and III' where only the singlet- 
likes states are mixed. It is therefore important to see how 
the probabilities of finding the system in the different states 
behave for long times after the evolution stage. 

First, from the results in SeclIV A2| we see that P,,/(tj) = 1, 



with PijiTs) the probability of finding the system after the 
evolution in the same state as the initial state, |<1>,-^. In gen- 
eral, Pij(Ts) decays to a saturation value in a time of the order 
ft/cru', this decaying time is sensitive to differences in the num- 
ber of '■'C atoms in each quantum dot leading to cr^ ^tcrj^. On 
the other hand the saturation value is robust because the limits 



3) = 1/2 and /|(tj 00) = 1/6 are independent 



of the variances as long as they are finite. For example, when 
\<^i) is singlet-like (excluding the trivial Pij(Ts) = 1 cases of 
Eq.(|63|) we get from Eq.( |62a[ i that P,;,(oo) - 1 /3. This means 
that the system is found with probability 2/3 in triplet-like 
states that will produce different outcomes after the joining 
stage. 

In summary, the evolution matrix is robust in two limits, 
at T5 — > 00 we get the saturation behavior while M^'^'(ts = 0) 
is just the identity matrix. For this reason we focus on f 
and Poo- In the following we discuss in detail the physics 
of the most interesting situations, those shown in FigjS] As 
discussed above, with the knowledge of k, j3 and any one of 
the LZ probabilities we can fully describe the separation and 
joining stages M^^' and M®. Whereas the disorder imposes 
the value of /3, the chosen detuning rates, v*J' and v*^\ con- 
trol both our reference probability (usually a LZ probability 
in a specific avoided crossing that we desire to maintain in the 
adiabatic limit) and the value of k. 

The most natural choice for /c is 1, when assuming that the 
detuning speeds at the separation and joining stages are the 
same. However that might not be the case experimentally. 
As an example, in Fig|9|b) we show that the simplest RC- 
charging/discharging circuit for pulse generation can induce a 
Ki^l situation. This happens for all cases if the boundary be- 
tween the (0,2) region and the (1,1) region [the detuning zone 
containing the avoided crossings, e 0, or more specifically 
fiiii > s > sii] is not centered between the preparation and 
measurement points (see FigjTfa)). Experimentally the point 
"p" is often closer to the boundary, as shown in the picture 
the separation stage is faster than the joining stage leading to 
a: < 1. This means that all the avoided crossings would be- 
come more effective at the joining stage and could, therefore, 
modify the outcome of the experiment. As we see below a sit- 
uation with Ki^l here has more impact on the final result than 
in 2DEG-based double dots where there is a unique energy 
gap due to interdot tunneling. 



1. Return probability when preparing tfie (0,2) ground state 

Here we study the return probability if the ground state 
is prepared, = l^l^!^'^. The goal is understanding how 
the return probability experiment is affected by the doublet- 
flipping interdot tunneling, i.e., for /3 > 0. We also investigate 
the effects of a mismatch in the rising times of the separa- 
tion and the joining stages: the behavior for different values 



of K. We start by choosing equal separation and joining de- 
tuning speeds, k-1. For any given value of /3 the speed cho- 
changes the outcome of the experiment. From 



sen, V 



,(s) . 



Fig|8[a), we see that the probability of state conversion at the 
avoided crossing II (that arises due to the disorder) is mean- 
ingful for the experiment, therefore, instead of showing the 
results as a function of v*J' and /3 we do it as a function of Pjj^ 
and p. The result presented in Fig|9|a) is the saturation return 
probability at the measurement region, f ™; after a joining 
stage that takes the system to a large measurement detuning, 
fim, so that £„ > fin- 

In order to study a region in the two dimensional space 
(p|j',y6) meaningful for the experiment we need an additional 
constrain: we will assume that at the time of performing the 
experiment one assures that the doublet conserving avoided 
crossing I, having the biggest gap in the system (for /5 < 1), 
is in the adiabatic limit. As seen in Fig|8ja) this guarantees 
that the separation stage is fully effective leaving the system 
m the (1,1) state, j^":"). The condition for adiabatic sepa- 
ration at I is equivalent to requiring a small complementary 



probability or leakage, P^^\ = f J ' » 0. From the LZ formula 



of Eq.([32]i and from the relations between the different gaps 
given in Eq.(|44d[), we obtain 



(64) 



In Figjoja) we present three curves Pf^ [p, Pf^^^ with I ap- 



proaching the adiabatic regime as P\^^]^ = 10 ,10 and 10 
reducing o 

detuning speed as 



It is clear that for reducing f j ., one needs to slow down the 



I old 



y 



leak 1 1 



With vo being the detuning speed that assures P\ 



,(s) 
leak 



10 



-2 



the cases with Pj^^j, = 10 and 10 shown in the figure are 
performed at speeds v'o/2 and vo/8, respectively. The smaller 
^leak' '■^^ more the experimental curve shifts to the left; such 
a shift increases the effectiveness of the avoided crossing II 
leading to a larger return probability. 

Having identified the experimentally relevant trajectories, 
we compute the probability of finding the system in a (0,2) 
state, for Tj = and for Tj — > 00, choosing P\^j^^ - 10" . 
The results, in general, only depend on [5 and k. As above 
described, a regular joining stage moves the detuning all the 
way to the > sn region, i.e., the measurement region; this 
leads us to the probabilities P™ and P". However, as sketched 
in Fig[8ja), one can also envision an experiment with a join- 
ing stage having a final detuning in the intermediate region 
given by ei < £,n < fiu and measuring the probabilities Pj, 
and P'oo- We obtain the intermediate final state simply by ex- 
cluding M^^iJ in the composition of M*j' given in Eq.(|56]i; the 



return probabilities are found applying Eq.p2|), i.e., summing 
the probabilities of ending in an (0,2) state. 

Here this distinction is very useful, because by studying 
how the P™ probabilities differ from the P' ones, we manage 
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FIG. 9. Return probability measured if the (0,2) ground state is prepared, see Fi^a). (a) Saturation return probability f'^ at the measurement 
region > en] as a function of fi and of the LZ probability in the DFAC II, Pjj', for k=\. The three trajectories fall in the experimental 

relevant regime: the LZ probability at I is assured to be in the adiabatic limit, i.e., p\ ' is small. The detuning velocity vo is chosen such that 

p[ ' = 0.01. As the speed decreases (vo/2 and Vo/8) 11 becomes more effective; both Pjj' and grow, (b) Example of Ki^ 1 for detuning pulses 
generated using a RC-charging/discharging circuit. If the zero detuning region is closer to the preparation point than to the evolution point 

(this is X < 1/2), a k^q < 1 situation is expected. Return probabilities along a physical relevant trajectory (-P}^*,^ = = 10"^) are show as a 
function of y6 and k in panels (c) and (d) [at the intermediate region, S\> era> £11], and (e) and (f) [at the measurement region]. The values for 
Tj = (panels (c) and (e)) are yS-independent and they remain close to 1 unless the joining stage becomes too rapid (/(■ = 9, see text). Panels (d) 
and (f) show the saturation return probabilities. In (d) the joining stage does not include the crossing II and P^^ =s 1/3 independently of j8 and k 
(for K < 2). In (f) the system passes through the crossing II before entering the measuring region and > 1/3. The return probability in (f) 
grows as the triplet-like states can return to (0,2) more effectively, i.e., the bigger it is /? and the slower it is the joining stage. 



to isolate the elfect of the avoided crossing II. Moreover, the 
return probabilities Pj, and (as shown in Fig|9|c) and (d)) 
behave as the well known spin-only double dot which gives 
Pj) = 1 and P'„ = 1 /3 i^^^'' This result is expected because, as 
we have shown above, the hyperfine field in the space of each 
Kramers doublet (for each dot) is isotropic in average, and 
also because the behavior of crossing I is completely analo- 
gous to the crossing involving the (1,1) and the (0,2) spin- 
singlets in 2DEGs' dots. The independence with /3 of the P' 
probabilities is trivial because for a given value of /5 at the 
separation stage the detuning speed is adjusted to reduce the 
leakage making the avoided crossing I adiabatic, and because 
the reduced joining stage does not involve the crossing II. 

Regarding the dependence with k of the probabilities P\ we 
find a drastic reduction when the joining stage is performed 
too fast, see for example k = 9. This means that the the avoided 
crossing I becomes non-adiabatic and the (1,1) singlet-like 
state is not effectively converted to an (0,2) singlet-like. By 
a reasoning similar to the one leading to Eq.(|65]l one gets. 

Because ^=;;JJ ^ = (^11)' ■ (66) 



Having P[^^^ = 10"^ we obtain that P® j, ^ 0.464 for k = 9, 
since the leaking at I becomes very high such too rapid join- 
ing stages are always avoided when designing the experiment. 
In what follows we concentrate the discussion on the smaller 
values of k. 

We note that for a- < 2 the P' values are virtually unaffected 
by changing the detuning speed of the joining stage (either de- 
creasing or increasing it). On the other hand, the P™ probabil- 
ity (see Fig|9|f)) is always affected by k when the double dot 
system have p > 0.001 . The smaller k the larger it becomes P™ 
increasing the difference with the standard case of P'^ a; 1/3. 
Those effects are due to the avoiding crossing II because it 
provides a new path /or returning to (0,2) to the (1,1) triplet- 
like states that are blocked at the avoided crossing I (states in 
Eq.(j4T]i taking d = d' = -). Of course, the enhancement of 
the return probability grows with fi even for k-\. Assuming 
that I is adiabatic also at the joining stage (which is a good 
approximation for a- < 2) we can write the return probability 
simply as, 

{ts) = fc ir.) f^ (T,) + 3 ft {ts) ff (r,) (67a) 
P™ (r,) = P- (r,) + ( 1 - P'(T,)) P® (67b) 
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The return probability f" (r,) is just Eq.( |62a| i, the well known 
probability for spin-only double dots at zero magnetic field 
of being found in a singlet-like state having started the evo- 
lution in the same singlet-like state. The enhancement term 
in f™ (ts) is proportional to f and to the probability of be- 
ing found after the evolution stage in a (1,1) triplet-like state, 
(1 - P' (js)). Experimentally it is possible to explore if the sys- 
tem responds to this description by designing a joining pulse 
with different v'J^ at the two avoided crossings. This can be 
done for instance by fixing v'J^ at I slow enough to maintain 
^leik small and by repeating the experiment for different val- 
ues of v*J' at the avoiding crossing 11. 

To summarize this part, when the prepared state is the 
ground state the return probability experiment produces a very 
similar output to the case of a spin-only double dot. A map- 
ping between the two situations in the absence of valley was 
done in our previous publication."*"* Here the hyperfine dy- 
namics remains isotropic in average but the joining stage is 
affected by a disorder induced avoided crossing (11) that al- 
lows for an enhancement of the return probability because it 
unblocks the triplet- like (1,1) states. 



tion of an (0,2) state, therefore, the dot R referred here must 
be identified with the dot which in the PRE is prepared with 
two electrons. 

First we study the saturation return probability at the mea- 
surement region, P™, fixing k=1. For a given double dot /? is 
fixed and the separation and joining stages can be performed 
slower or faster changing the whole set of LZ-probabilities. 



With the assumptions discussed in Sec IV A 1 all the LZ- 
probabilities are linked once /3 and k are given. In this case, 
we choose the LZ probability and (3 as the two variables 
that define P™. We do so because here, as shown in Fig[8jb), 
the probability of state conversion at the avoided crossing IIP 



plays a crucial role. The result is shown in Fig 10 a) and (b) in 
the whole parameter space. In order to study situations exper- 
imentally relevant we assume that the detuning speed is tuned 
so that the avoided crossings with the biggest gaps are the adi- 



abatic limit, i.e., P, 



leak ' 



0. From Eq.&2\ and Eq.(|44d 



we find that the LZ probability of state conversion at III' is. 



^III'V^' leak/ ^ yie-dk) 



(68) 



In Fig 10 'a) and (e) we show two curves with the values of 
P\f,, resulting from choosing P{^' = 10"^ or, as it results from 



2. Return probability when preparing (0,2) higher energy states halving the speed of the separation Stage, 10 



We fix F 



10- 



The purpose of this section is to demonstrate the complica- 
tions that arise when the prepared state is not the (0,2) ground 



state. We show that the simple procedure followed in Sec IV B 
of reducing the leakage probability in the avoided crossings 
with the biggest gaps (for /3 < I the doublet conserving, I and 
I') does not assure a successful separation stage. After a suc- 
cessful separation stage the system is in a (1,1) configuration, 
therefore the probability of finding the system in (0,2) at the 
evolution region, P*^, must be close to zero. 

The situation we study is sketched in Figj8|b), the prepared 
state is the highest energy excited state, \^'p^ - \sff'^y As 
opposite to the ground state case, the separation stage first in- 
volves the LZ process in III', which is an avoided crossing 
induced by disorder (a DFAC). From Eq.(|44]), we have that 
^iir = "V^^is this links the probabilities of state conversion in 
r and III' (assuming the detuning speed is the same) as shown 
in Fig|5]^b). For weak disorder Am' is small and therefore the 
adiabatic condition cannot be met while still having the pas- 
sage time larger than the characteristic time of the HFI. In- 
stead, we explore the effect of III' being a diabatic process in 
general, although the limit of III' being adiabatic is explored 
when yS is close to 1 and Anr ~ Ay . 

Once an excited state is prepared multiple Landau-Zener 
processes — including LZ loops as discussed in Sec IV A 1 — 
affect the separation and joining stages. Furthermore, the 
Landau-Zener sequence of relevant avoided crossings differs 
for A*^ < A^ and A*' > A*^. Obviously, these two cases can be in- 
vestigated experimentally in the same device first by perform- 
ing the PRE after the preparation of a desired excited state 
in the (0,2) configuration and then repeating the measurement 
preparing an equivalent excited state in the (2,0) charge con- 
figuration. In our calculations we always assume the prepara- 



leak ~ ' ™ intermediate situation between 
the latter two conditions) and we investigate the effectiveness 
of the separation sta ge b y plotting as a function of yS; pan- 
els (c) and (f) in Fig 10 correspond to A*^ < A*^ and A^ > A^, 
respectively. We can see that for (3 < 0.001 the separation is 
effective, P^ » 0. The system is essentially clean and only the 
DCAC r is active. This is consistent with the result shown 
in Fig|5];b), for very small yS the LZ process in the DFAC can 
be neglected even though the LZ processes in the DCAC (vis- 
ited with the same detuning speed) approaches the adiabatic 
regime. This can also be seen from Eq.(|68]l, the small value 
of Pj^l^j^ — curves presented in Fig 10 a) and (e) — implies only 
for j6 very small that P[jj, can be neglected. 

As yS grows up to 0.1 P^ also grows: the separation stage 
becomes ineffective. This results from a non-negligible Pjj*, . 
However, for the range of fi between 0.1 and 1, P*^ decreases 
for A*^ < A^ and the separation stage improves, whereas for 
^ pe i^ggpg growing with yS. This is caused by the ar- 
rangements of the LZ processes, as sketched with grey arrows 
in the subpanels (i) and (iii) of Fig|7|c). A finite PjjJ, implies 
the passage through the DCAC IV' m a detuning/energy con- 
figuration that differs depending on A** and A*". Due to the 
gap ratios in Eq.(|44]) the LZ processes at IV' and IV are al- 
most adiabatic (from Pj^^,^ = 10"^ we get Pj^, ~ 0.969). This 
means that for A*^ < A*' the state I^L''^''^ is converted at IV' 
into \S 'l'*!^*^ and IV is irrelevant, whereas for A*' < A*^ both IV' 

and IV are involved leading to the state 

In the limiting case of p - 1, the LZ processes in all the 



,(s) 



1 



avoided crossings can be considered adiabatic. Since Pj 
the system follow 100% the path indicated by the grey arrows 
in subpanels (i) and (iii) of Fig|7]^c). This means that for A*' < 
A^ the system leaves IV at state 15*^'!''^ and then fully follows 
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(see Fi^b) and Figj6jfor the evolution and measurement regions, and respectively). 
Panels (a), (b), (c) and (d) are for < A"^ while panels (e), (f), (g) and (h) are for A'' > A'*; see Fig|7|c) and Eq.j56[(. Panels (a) and (e) show 
the saturation return probability at the measurement region as a function of [} and of the LZ probability in the disorder induced crossing 
Iir, PjjJ, , for = 1. The two trajectories indicate the experimental conditions assuring that the LZ probability at crossing F approaches the 



adiabatic limit, i.e., P^, 



is small. We choose P\,' = 10"' and present the results as a function of /? for different values of k in the six panels (b) 
and (f), [probabilities at region R^, P^], (c) and (g) [probabilities at region R™ for r, = 0, P™], and (d) and (h) [probabilities at region R™ for 
Ts — > 00, P™]. First, P^ is not zero and it depends on y6, i.e., the system can remain in (1,1) after separation stage. Moreover, the ineffective 
separation stage combined with an ineffective joining stage make the return probabilities for = smaller than 1 . In panel (d) we find that the 
saturation return probability can be around 1/6 for a broad region in parameter space (/?, k), B; the associated r, = value in panel (c) is around 
1/2. Results in (g) and (h) for A'' > A"* do not show the 1/6 saturation value with the same robustness as a function of /? and a-; such a difference 
could be tested by repeating the experiment preparing an (2,0) state (see text). 



Ill (see Fig|6j) ending up in the \S "i"!^^^ (0,2) state. This justifies 
the fact that for /? = 1 the separation stage is fully ineffective 
for < and P*^ — > 1 . On the other hand, the separation 
stage is effective, P*^ — > 0, in the adiabatic limit for A"^ < A^: 
the system leaves IV' at state jS'^'l^*^ and after the adiabatic 
state conversion at IF (see Figj6j) it ends up in the j^'^'l''^ state. 

It is clear, from Fig|6j that only when the ground state is 
prepared the first avoided crossing involved in the separation 
stage is a doublet conserving one. By the present example, 
we have shown a feature common to the case of preparing any 
excited state of (0,2), namely, that if the region of detuning 
involving the avoided crossings is visited with a constant de- 
tuning speed the separation stage is, in general, not effective. 
A return probability experiment performed in such a condition 
does not achieve the first conceptual goal of the RPE: separat- 
ing the electrons. A detuning pulse engineered with different 
speeds at the different crossings can be the solution. This is 
the most important conclusion of this section. 

For the sake of completeness, despite the ineffectiveness of 
the separation stage, we now show what return probabilities, 
P™, would be observed at the measurement region as a func- 



tion of p and fixing Pf^^^^ - 10"^. The results are shown for 
different values of k in panels (c), (d), (g), and (h) of Fig 10 
Similarly to the case in the previous section, k - 9 implies 
a too fast joining stage that affects even the clean limit. For 
the smaller values of k the first observation is the trivial re- 
cover foryS ^ of the P™ 1, P™ ss 1/3 — consistent with the 
results with no valley mixing.'*'* In a very disordered system 
with y6 — > 1 and A"* < A^, as discussed above, the separation 
stage is effective and again we recover the P™ » 1, P™ a; 1/3 
physics. On the other hand, for yS — > 1 and A'' < A^ one 
measures P™ ^ 1 and P™ ^ 1, but the system does not have 
electrons separated in (1,1) states in the evolution stage. 

Notably, for A"* < A'' we find an intermediate region in the 
P axis, see zone B in Fig 10 d), for which P™ ^ 1/6; the result 
persists even when the velocities of the separation and joining 
stages differ (for k < 2). This value agrees very well with the 
saturation probability, 0.17, measured in Ref. 14 However, 
whereas in the experiment for t, — > the return probability 
is 1, in our case study the ineffective separation and joining 
stages produce, in region B, P™ » 1/2 as shown in Fig 10 c). 

Here, the robustness of the P™ » 1 /2, P™ » 1/6 physics can 
be understood because the non-adiabatic LZ probabilities are 
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linked by, f*? = 1 - (1 - with h = s,j. We find that 

the region B approximately cover the values of /3 that produce 
0.3 > Pjjj, > 0.8. We do not discuss further this situation be- 
cause, as mentioned above, it is not a proper return probability 
experiment. However, a final interesting remark can be made 
regarding the experimental results of Ref.[14l If the measured 
P™ ~ 1/6 were due to the preparation of the highest excited 
state (0,2), as described here, then the preparation of the (2,0) 
analog state would not show such a saturation value with the 
same robustness as a function of k and this is shown in 
FigfTolh) withAL<A'*. 



C. RPE in the high magnetic field limit 

In spin-only double quantum dots one of the most com- 
mon procedures is applying an external magnetic field much 
larger than the hyperfine interaction energy. Focusing on the 
(1,1) states, the Zeeman interaction shifts in energy the spin 
polarized triplets leaving in a double degenerated subspace 
the spin-singlet state \S} = |Ti)-|iT) and the triplet-0 state 
\To} - |Ti) + liT)- This two level system defines a qubit that 
can be manipulated. The return probability experiment be- 
comes simpler than for the case for zero magnetic field. The 
(0,2) singlet is prepared and after the separation stage the in- 
homogeneous part of the hyperfine interaction mixes the (1,1) 
singlet with the triplet-0. The result of averaging over many 
realizations leads to P^o = 1 /2 in the high detuning limit (neg- 
ligible tunneling exchange) .I^^MD 

The search for two-level subspaces that can be manipulated 
naturally leads to the investigation of the effect of external 
magnetic field. Here we focus on how the return probabil- 
ity behaves in such cases for double dots in nanotubes with 
valley mixing. The finite field case was studied for the clean 
nanotube system in Ref. 44; Here, our goal is exploring the 
high magnetic field case, i.e., when the magnetic energies are 
much bigger than the hyperfine energy, the spin-orbit splitting 
and the valley mixing energies. The form of the single-particle 
solutions at those limits determines the nature of the (0,2) and 
(1,1) solutions and the behavior of the experiment. 

We solve the dot Hamiltonian of Eq.Q, in Fig{TT|a) and 
(b) we plot the single-particle energies as a function of the 
Zeeman energy for the isolated left and right dot assuming 
they have different valley mixing energies. We present in the 
same panels the cases with B || z and B ± z where z is taken 
along the tube axis. The two particle energies for the (0,2) and 
the (1,1) charge configurations (at f = or in high detuning 
limit) are plotted in panels (c) and (d) of FigfTT] they follow 
from Eq.([T7]i and Eq.([T8]). 

We first focus on the B || z case, neglecting the HFI the 
spin projection along the z-axis, cr, is a good quantum num- 
ber. The general single-dot single-particle solutions for this 
case is given in Eq.(|5]l. The solutions in valley space are the 

spinors j^"^, they point parallel {+) and antiparallel (-) to 

the valley vector g^. In the high-field limit the valley vector 
component out-of-the-plane dominates due to the diamagnetic 
component that couples with the valley degree of freedom. 



|5o-| = IcrAfo - Zioi-bl » a|.^,, and the valley spinors ^g^., +^ be- 
come |t), with T = /r, K' (the ±1 eigenstates of t^). This means 
that the solutions for the dots are simply If rcr) with eigenen- 
ergies. 



E^ra- - (E^a- + EorbT - A,so TO") /2, 



(69) 



where r is -i-l or -1 for K or K', respectively. This result 
is identical to a clean nanotube quantum dot with spin-orbit 
coupling in a parallel field. This is confirmed by the (0,2) 
and (1,1) spectrum at the high B || z limit in Fig[TT|c) and (d). 
The available return probability situations correspond to those 
reported in Ref 44, leading to P^o = 1/2 or P^o = 1 depending 
on the prepared state. No disorder-induced avoided crossings 
survive in such a limit because the solutions in both dots have 
the same spin and valley characteristics. Any given (0,2) state 
is mixed, due to interdot tunneling, with a single (1,1) state 
and the gap energy is 2 ^/2\t\ in all six avoided crossings. 

We now focus on the perpendicular magnetic field case, we 
take B = B^ii. We choose the single particle basis l^Ts^} with 
■s.v =Ta J Ix (or if = +) the spin projection along the direction of 
the magnetic field. The spin-orbit coupling term in the Hamil- 
tonian of Eq.([T]i is -^AfoTscr. and mixes l^Ks^} with \^Ks^}. 
Such a mixing can be neglected in the high magnetic field 
limit because states with opposite spin projections are split 
in energy by the Zeeman term, EsTqct^. However, the valley 
mixing term in the Hamiltonian mixes states of opposite val- 
ley but identical spin projection s^', the effective Hamiltonian 
at high Bjt for the electrons s ^ in dot f is the valley operator, 



E cS y 



4jf'exp(-i^l^,) 
The solutions are the in-plane valley spinors 



(70) 



(71) 



with = ±. So the full single particle solutions are |f m^,? 
If) ® l^f) ® k.i), with energies. 



(72) 



The four (1,1) lowest energy states FigjTTJd) are the (1,1) 

four Slater determinants with =i.i, i.e., with mr = 

± and Ml = +. In general, A^^, A^^, and therefore the 
valley mixing breaks the degeneracy between these four states 
because (disregarding a global energy shift) the energies are. 
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KK' 



(73) 



This is a clear distinction of this system with spin-only double 
dots, there is not any double degenerated subspace (analogous 
to the one spanned by ITq) and \S)) in the (1,1) configuration. 
The last issue we investigate for perpendicular field is how 

the (0,2) ground state, mixes with the (1,1) states. This 
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FIG. 11. Qualitative differences between parallel and perpendicular (to the tube axis) magnetic fields and their interplay with disorder. We 
take a common spin-orbit splitting Aj^, the disorder generates dot-dependent valley mixing energies, Aj^j^., 4^ A^j^, . Panels (a) and (b) show 
the single-particle eigenvalues for the isolated right and the left dot, respectively, as a function of the Zeeman energy we distinguish the 
parallel from the perpendicular magnetic field case. At the high-field limit, A^o |-£s|, the valley mixing effect is irrelevant for the parallel 
case (due to the diamagnetic term) while on the other hand it is important for the perpendicular case. It follows that in the high-field limit 
parallel magnetic fields leads to the same physics than in a clean nanotube.'*'' In panels (c) and (d) we plot the energies of the (0,2) and (1,1) 
eigenstates (in the high detuning limit or for r = 0). Differences in the two valley mixing energies generate splitting of the four lowest energy 
(1,1) states for perpendicular magnetic field case. In panel (e) we plot the mixing between (0,2) and (1,1) states as a function of the detuning 
for a case in which the Zeeman energy due to the perpendicular magnetic field is dominant. We focus on the mixing involving the (0,2) ground 
state, GS(0,2). The (0,2) ground state is spin polarized and therefore can only mix with the four lowest energy states of (1,1), L-, R-, that are 
also spin polarized. As follows from Eq.|74| the tunneling is determined by the valley properties and therefore the phase difference between 
the valley mixings, ifg^i^, tunes the mixing. In the insets above we show that, unless i/j/^r = rm, one should expect mixing of GS(0,2) with the 
four L-, R- states of the (1,1) configuration. 

Therefore, for = 0, sin ^ = and the (0, 2) ground state 
mixes only with the two intermediate energy states, while on 
the other hand for ipjn^ - n, cos ^ - and the mixing is only 
with the other two states. For other values of tpjn there are 
tunneling amplitudes with the four states. 

All these effects have to be considered if a return probabil- 
ity experiment is to be performed preparing the (0,2) ground 
state. We note that if ^ |A^^, - A^^,| » cth, the hyperfine 
interaction would not induce mixing between the (1,1) states. 
However, the experiment, if performed with different detun- 
ing speeds and different waiting times Tj, can provide infor- 
mation about the values of A|.^, and (fRi. Some oscillations 
in the signal would not average out when repeating the single- 
shot cycle if the values of A^^, are fixed in time. This would 
be the case if its origin is indeed static disorder and defects. 

V. CONCLUSIONS 

The return probability experiment is a standard tool for 
characterizing dephasing in double quantum dots systems. 
The valley degree of freedom in nanotubes introduces more 
states in the system. The spin-orbit coupling, the Zeeman in- 
teraction, and the diamagnetic effect of a parallel magnetic 



is important for the return probability experiment if the ground 
state is to be prepared. First due to spin conservation of the 
single particle tunneling, Ht = -I^iTqctq, the two spin-down 
electrons can only mix with the four lowest energy (1,1) states 
having the same spin configuration. The result is confirmed in 
FigjTTJe) where we plot the two particle spectrum as a func- 
tion of detuning choosing a high value of B^. We notice that 
the phase difference between the valley mixing interactions in 
the two dots, -<f^f^, -v'kk'^ plays a key role in weighting 
the tunneling between the (0,2) ground state and the different 
(1,1) states. This is shown in the insets in Fig[TT|e). 

Such a dependence follows from computing the matrix ele- 
ments of the tunneling Hamiltonian between solutions with a 
given spin in different dots, we get: 

/T I 17 ID \ f-fCOS^ forMR = ML = + 

I -ft sin forMR = -ML = ± 

Then the application of the tunneling Hamiltonian to the (0,2) 
ground state gives , 
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field, are elements to be considered when performing the ex- 
periment. In ultra clean nanotube quantum dots, in which one 
can neglect valley mixing effects, the return probability exper- 
iment presents many different scenarios. In such a case, the 
origin of the different cases lies in the effective hyperfine in- 
duced dynamics within subspaces of (1,1) states. Those sub- 
spaces are defined by the chosen parameters. On the other 
hand, the separation and joining stages of the experiment are 
trivial. 

In this paper we focus on nanotubes for which the valley 
mixing energy is much larger than the hyperfine energy, cru, 
and therefore it cannot be neglected. The presence of valley 
mixing could be due to defects and impurities in the nanotube. 
Since the disorder profile is position dependent the resulting 
valley mixing terms are dot-dependent. First, we focused on 
the zero-magnetic field case and solved the single-dot prob- 
lem. At the dot ^ = L, R the spin-orbit coupling and the valley 
mixing split the four states, |^d^cr^, an energy d^A^/2 into two 
Kramers doublets, the low energy one, d^ = - and the high en- 
ergy one df = -. Because cru ^ the hyperfine interaction is 
unable to mix states from different doublets. We have shown 
that the effective hyperfine interaction in the space of the two 
states in a Kramers doublet remains statistically isotropic as 
for the case of clean nanotubes. The dynamics due to hy- 
perfine interaction in the (1,1) states is, as in spin only double 
dots, dictated by precession along the local hyperfine effective 
fields that are different in the two dots and in the two Kramers 
doublets. 

In the high detuning limit in absence of hyperfine interac- 
tion the sixteen (1,1) eigenenergies split into the four fourfold 
degenerated subsets, (LdL, RdR). The action of the hyperfine 
interaction cannot change the Kramers doublet index d^ and 
the final state lies in the same subset that the initial state. We 
get the standard behavior for the hyperfine dynamics within 
each subset: if starting in a singlet-like state, the probability 
to be found at » S/ctjj in the same state saturates (aver- 
aging over many single-shot experiments) to a value 1/3. It 
is worth noting that the isotropic hyperfine field we have ob- 
tained is the most effective coherent form for dephasing the 
system. This isotropy, in average, leads to the lowest possible 
saturation value; the presence of anisotropy in the hyperfine 
field — that can arise in clean systems for perpendicular mag- 
netic fields — increases the saturation value up to 3/8."*"* The 
situation we are studying is therefore the most favorable to 
achieve a secondary goal, obtaining the lowest possible satu- 
ration values allowed in the system given that the only exper- 
imental result available is lower than 1 /3.G3 

Whereas the valley mixing does not introduce qualitative 
changes in the hyperfine interaction it does so for the single- 
particle interdot tunneling. The solutions in each dot point in 
noncolinear directions in valley space. This is shown to pro- 
duce nonzero tunneling matrix elements between solutions in 
different dots with the same spin cr irrespective of the Kramers 
doublet index. In the two particle spectrum reflecting the mix- 
ing between the (1,1) and (0,2) states as a function of de- 
tuning, we show that four new avoided crossings appear as 
a consequence of the doublet-flipping interdot tunneling am- 
plitudes. These avoided crossings can play a crucial role in 



the separation and joining stages modifying the output of the 
experiment. 

For gaining qualitative knowledge of the implications of the 
novel crossings we adopt a low tunneling picture that allows 
us to use the Landau-Zener formula for each avoiding cross- 
ing. In such an approximation the ratio between tunneling 
amplitudes squared, /3, defines the relations between all the 
gaps in the avoided crossings. We show that when any ex- 
cited (0,2) state is prepared the separation stage first visit an 
avoided crossing induced by disorder. Any separation stage 
that crosses the tunneling mixing region with the same detun- 
ing speed (a not fine tuned detuning pulse) can lead to im- 
perfect separation leaving the electrons in the same dot. This 
is not consistent with a proper return probability experiment 
since it leads to f (tj = 0) < 1 ; for example we show that when 
preparing the highest exited state one gets P{ts - Q) ~ Ijl 
and Poo ~ 1/6 in a broad region in parameter space. Despite 
the agreement of the saturation value with the experimental 
one,™ 0.17, there exist a strong discrepancy at short evolv- 
ing times because the measurement goes to 1 . As mentioned 
above more experimental work would be desirable, including 
the dependence with magnetic fields. 

We also showed that preparing the ground state leads to a 
separation stage that first visits the doublet conserving avoided 
crossing 1. This crossing can be visited adiabatically and no 
further avoided crossings are involved in the separation stage. 
We find that the return probability experiment behaves very 
similar to the spin-only case, except that the (1,1) triplet-like 
states (with 2/3 probabilities of being occupied in average at 
Tj — > oo) have a nonzero probability, P®, of returning to (0,2) 
due to the disorder induced LZ-process at the avoided cross- 
ing 11. For the saturation return probability this implies that it 
cannot be smaller than 1/3, namely, x 1 /3-(-2/3P®. 

Finally, we have studied the high magnetic field limit. 
When the magnetic field is parallel to the tube axis the dia- 
magnetic effect competes directly with the valley mixing term 
because both operate in valley space for electrons with the 
same spin projection. This makes the valley mixing irrelevant 
for high magnetic fields, the single particle solutions become 
valley polarized and the return probability experiment pro- 
duces the same results predicted for clean nanotubes in paral- 
lel magnetic fields."*"* On the other hand, if the magnetic field 
is perpendicular to the tube axis there is no diamagnetic inter- 
action and the valley mixing splits states with the same spin 
that otherwise would be degenerated at high magnetic fields. 
We show that the phase difference between the valley mixing 
terms in the two dots controls how the interdot tunneling be- 
haves and this changes the avoided crossing pattern between 
the (0,2) ground state and the (1,1) four lowest energy states 
(all being fully spin-polarized). Though the hyperfine inter- 
action would play no role, variations of the return probability 
experiment could give insight of the valley mixing terms by 
studying the oscillations in the signal. The combination of the 
multiple mixing of the (0,2) ground state with the fact that no 
double degenerate (1,1) subspace remains due to disorder — 
analogous to the singlet and triplet-0 subspace in spin-only 
double dots — can be important in the interpretation of Pauli 
blockade measurements in nanotubes in a perpendicular mag- 
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